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\& • This overview focuses on the notion of partial dynamical symmetry (PDS), for which a pre- 

scribed symmetry is obeyed by a subset of solvable eigenstates, but is not shared by the Hamil- 
tonian. General algorithms are presented to identify interactions, of a given order, with such 
intermediate-symmetry structure. Explicit bosonic and fermionic Hamiltonians with PDS are con- 
structed in the framework of models based on spectrum generating algebras. PDSs of various 
types are shown to be relevant to nuclear spectroscopy, quantum phase transitions and systems 
with mixed chaotic and regular dynamics. 
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1 Introduction 



Symmetries play an important role in dynamical systems. Constants of motion associated with a 
symmetry govern the integrability of a given classical system. At the quantum level, symmetries provide 
quantum numbers for the classification of states, determine spectral degeneracies and selection rules, 
and facilitate the calculation of matrix elements. An exact symmetry occurs when the Hamiltonian 
of the system commutes with all the generators of the symmetry-group G, [H , g t ] = 0. In this 
case, all states have good symmetry and are labeled by the irreducible representations (irreps) of G. 
The Hamiltonian admits a block structure so that inequivalent irreps do not mix and all eigenstates in 
the same irrep are degenerate. In a dynamical symmetry the Hamiltonian commutes with the Casimir 
operator of G, [ H , Gq ] = 0, the block structure of H is retained, the states preserve the good symmetry 
but, in general, are no longer degenerate. When the symmetry is completely broken then [H , g^] ^ 0, 
and none of the states have good symmetry. In-between these limiting cases there may exist intermediate 
symmetry structures, called partial (dynamical) symmetries, for which the symmetry is neither exact 
nor completely broken. This novel concept of symmetry and its implications for dynamical systems, in 
particular nuclei, are the focus of the present review. 

Models based on spectrum generating algebras form a convenient framework to examine underly- 
ing symmetries in many-body systems and have been used extensively in diverse areas of physics [I]. 
Notable examples in nuclear physics are Wigner's spin-isospin SU(4) supermultiplets SU(2) single- 
j pairing [3J, Elliott's SU(3) model [I], symplectic model [5], pseudo SU(3) model [5J, Ginocchio's 
monopole and quadrupole pairing models [7J, interacting boson models (IBM) for even-even nuclei [8] 
and boson- fermion models (IBFM) for odd-mass nuclei [9]. Similar algebraic techniques have proven 
to be useful in the structure of molecules and of hadrons [12]. In such models the Hamiltonian 

is expanded in elements of a Lie algebra, (Go), called the spectrum generating algebra. A dynamical 
symmetry occurs if the Hamiltonian can be written in terms of the Casimir operators of a chain of 
nested algebras, Go D Gi D . . . D G n [13] . The following properties are then observed, (i) All states 
are solvable and analytic expressions are available for energies and other observables. (ii) All states are 
classified by quantum numbers, \a$, a\, . . . , a n ), which are the labels of the irreps of the algebras in the 
chain, (iii) The structure of wave functions is completely dictated by symmetry and is independent of 
the Hamiltonian's parameters. 

A dynamical symmetry provides clarifying insights into complex dynamics and its merits are self- 
evident. However, in most applications to realistic systems, the predictions of an exact dynamical 
symmetry are rarely fulfilled and one is compelled to break it. The breaking of the symmetry is 
required for a number of reasons. First, one often finds that the assumed symmetry is not obeyed 
uniformly, i.e., is fulfilled by only some of the states but not by others. Certain degeneracies implied 
by the assumed symmetry are not always realized, (e.g., axially deformed nuclei rarely fulfill the IBM 
SU(3) requirement of degenerate /3 and 7 bands [5]). Secondly, forcing the Hamiltonian to be invariant 
under a symmetry group may impose constraints which are too severe and incompatible with well- 
known features of nuclear dynamics (e.g., the models of [7J require degenerate single-nucleon energies). 
Thirdly, in describing transitional nuclei in-between two different structural phases, e.g., spherical and 
deformed, the Hamiltonian by necessity mixes terms with different symmetry character. In the models 
mentioned above, the required symmetry breaking is achieved by including in the Hamiltonian terms 
associated with (two or more) different sub-algebra chains of the parent spectrum generating algebra. 
In general, under such circumstances, solvability is lost, there are no remaining non-trivial conserved 
quantum numbers and all eigenstates are expected to be mixed. A partial dynamical symmetry (PDS) 
corresponds to a particular symmetry breaking for which some (but not all) of the virtues of a dynamical 
symmetry are retained. The essential idea is to relax the stringent conditions of complete solvability so 
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exact symmetry degeneracy 
all states with good symmetry 

[H, gi } = 9l eG 



dynamical symmetry splitting no mixing 



all states solvable, with good symmetry 
Gq D Gi D . . . D G n 

|ao, ai, . . . , On) 

[H,C Gk }=0 H = Ea k C Gk 



partial symmetry 



part of the degenerate states (solvable) with good symmetry 

[H, gi ] f o 
[H, m m = 9l eG 



partial dynamical symmetry 



part of non-degenerate states (solvable) with good symmetry 

[H,C G ] f 

[h, c G m = 



broken symmetry 



none of the states solvable, without good symmetry 
Figure 1: Hierarchy of symmetries. 



that the properties (i)-(iii) are only partially satisfied. It is then possible to identify the following types 
of partial dynamical symmetries 

• PDS type I: part of the states have all the dynamical symmetry 

• PDS type II: all the states have part of the dynamical symmetry 

• PDS type III: part of the states have part of the dynamical symmetry. 

In PDS of type I, only part of the eigenspectrum is analytically solvable and retains all the dynamical 
symmetry (DS) quantum numbers. In PDS of type II, the entire eigenspectrum retains some of the DS 
quantum numbers. PDS of type III has a hybrid character, in the sense that some (solvable) eigenstates 
keep some of the quantum numbers. 

The notion of partial dynamical symmetry generalizes the concepts of exact and dynamical symme- 
tries. In making the transition from an exact to a dynamical symmetry, states which are degenerate 
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in the former scheme are split but not mixed in the latter, and the block structure of the Hamiltonian 
is retained. Proceeding further to partial symmetry, some blocks or selected states in a block remain 
pure, while other states mix and lose the symmetry character. A partial dynamical symmetry lifts 
the remaining degeneracies, but preserves the symmetry-purity of the selected states. The hierarchy of 
broken symmetries is depicted in Fig. 1. 

The existence of Hamiltonians with partial symmetry or partial dynamical symmetry is by no 
means obvious. An Hamiltonian with the above property is not invariant under the group G nor does it 
commute with the Casimir invariants of G, so that various irreps are in general mixed in its eigenstates. 
However, it posses a subset of solvable states, denoted by in Fig. 1, which respect the symmetry. 
The commutator [H , gi] or [H , Cq] vanishes only when it acts on these 'special' states with good 
G-symmetry. 

In this review, we survey the various types of partial dynamical symmetries (PDS) and discuss 
algorithms for their realization in bosonic and fermionic systems. Hamiltonians with PDS are explicitly 
constructed, including higher-order terms. We present empirical examples of the PDS notion and 
demonstrate its relevance to nuclear spectroscopy, to quantum phase transitions and to mixed systems 
with coexisting regularity and chaos. 

1.1 The interacting boson model 

In order to illustrate the various notions of symmetries and consider their implications, it is beneficial 
to have a framework that has a rich algebraic structure and allows tractable yet detailed calculations 
of observables. Such a framework is provided by the interacting boson model (IBM) [5| ITU - TIE] , widely 
used in the description of low-lying quadrupole collective states in nuclei. The degrees of freedom of 
the model are one monopole boson (V) and five quadrupole bosons (dt). The bilinear combinations 
{s^s, s^d^, d \ x s, d^dfj,/} span a U(6) algebra. These generators can be transcribed in spherical tensor 
form as 

h s = s j s , ul L) = (Sd)f ] L = 0,1, 2, 3, 4 

n{?> = d\s + 8 % , IT?) = <(dt s - 8 %) , (1) 

where d M = (— l^cL.^, and standard notation of angular momentum coupling is used. U(6) serves 
as the spectrum generating algebra and the invariant (symmetry) algebra is 0(3), with generators 
1$ = y/WUjp. The IBM Hamiltonian is expanded in terms of the operators ([I]) and consists of 
Hermitian, rotational-scalar interactions which conserve the total number of s- and d- bosons, N = 
n s + ha = s^s + ^ d^dfj,. Microscopic interpretation of the model suggests that for a given even-even 
nucleus the total number of bosons, N, is fixed and is taken as the sum of valence neutron and proton 
particle and hole pairs counted from the nearest closed shell [TT] . 
The three dynamical symmetries of the IBM are 

U(6) D U(5) D 0(5) D 0(3) anharmonic spherical vibrator 

U(6) D SU(3) D 0(3) axially— deformed rotovibrator (2) 

U(6) D 0(6) D 0(5) D 0(3) 7— unstable deformed rotovibrator 

The associated analytic solutions resemble known limits of the geometric model of nuclei [18], as in- 
dicated above. Each chain provides a complete basis, classified by the irreps of the corresponding 
algebras, which can be used for a numerical diagonalization of the Hamiltonian in the general case. 
In the Appendix we collect the relevant information concerning the generators and Casimir operators 



5 



of the algebras in Eq. (J2J. Electromagnetic moments and rates can be calculated in the IBM with 
transition operators of appropriate rank. For example, the most general one-body E2 operator reads 



T(E2) = e B [U^+ X U^] 



(3) 



A geometric visualization of the model is obtained by an energy surface 



(p,r,N\H\f3, T ,N) 



(4) 



defined by the expectation value of the Hamiltonian in the coherent (intrinsic) state [T9jl20] 



\P,r,N) 



(N\r 1/2 (bl) N \0) , 

(1 + f3 2 )- l/2 [f3 cos 7 4 + P sin 7 (4 + dl 2 )/V2 + s f ] . 



(5a) 
(5b) 



Here (/3, 7) are quadrupole shape parameters whose values, (/?o, 70), at the global minimum of En(/3, 7) 
define the equilibrium shape for a given Hamiltonian. The shape can be spherical ((3 = 0) or deformed 
(f3 > 0) with 7 = (prolate), 7 = n/3 (oblate), 7-independent, or triaxial (0 < 7 < vr/3). The latter 
shape requires terms of order higher than two-body in the boson Hamiltonian [2TH22] . The equilib- 
rium deformations associated with the Casimir operators of the leading subalgebras in the dynamical 
symmetry chains (J2J) are, /3q = for U(5), (/3q = v2>7o = 0) for SU(3) and (/3q = 1, 70 arbitrary) for 



PDS of type I corresponds to a situation for which the defining properties of a dynamical symmetry 
(DS), namely, solvability, good quantum numbers, and symmetry-dictated structure are fulfilled exactly, 
but by only a subset of states. An algorithm for constructing Hamiltonians with this property has been 
developed in [23] and further elaborated in [21]. The analysis starts from the chain of nested algebras 



where, below each algebra, its associated labels of irreps are given. Eq. fl6]) implies that Gdyn is the 
dynamical (spectrum generating) algebra of the system such that operators of all physical observables 
can be written in terms of its generators [H1[T3]; a single irrep of Gdyn contains all states of relevance 
in the problem. In contrast, G S ym is the symmetry algebra and a single of its irreps contains states that 
are degenerate in energy. A frequently encountered example is G sy m = 0(3), the algebra of rotations 
in 3 dimensions, with its associated quantum number of total angular momentum L. Other examples 
of conserved quantum numbers can be the total spin S in atoms or total isospin T in atomic nuclei. 

The classification ((6]) is generally valid and does not require conservation of particle number. Al- 
though the extension from DS to PDS can be formulated under such general conditions, let us for 
simplicity assume in the following that particle number is conserved. All states, and hence the rep- 
resentation [h], can then be assigned a definite particle number N. For N identical particles the 
representation [h] of the dynamical algebra Gdyn is either symmetric [N] (bosons) or antisymmetric 
[1^] (fermions) and will be denoted, in both CclSGS, clS 

[Tijv]- For particles that are non-identical under 
a given dynamical algebra Gdyn, a larger algebra can be chosen such that they become identical under 



0(6). 



2 PDS type I 




d a 



A 



sym 



(6) 
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this larger algebra (generalized Pauli principle). The occurrence of a DS of the type signifies that 
the Hamiltonian is written in terms of the Casimir operators of the algebras in the chain 



H 



DS 



(7) 



G 



and its eigenstates can be labeled as |[^at](E) ... A); additional labels (indicated by ... ) are suppressed 
in the following. Likewise, operators can be classified according to their tensor character under ([6]) as 

T[h n ]{a}\- 

Of specific interest in the construction of a PDS associated with the reduction ([6]), are the n-particle 
annihilation operators T which satisfy the property 

f [hn]m \[h N )(V o )A) = 0, (8) 

for all possible values of A contained in a given irrep (So) of G. Any ra-body, number-conserving normal- 
ordered interaction written in terms of these annihilation operators and their Hermitian conjugates 
(which transform as the corresponding conjugate irreps) 



a,/3 



(9) 



has a partial G-symmetry. This comes about since for arbitrary coefficients, A a p, H' is not a G-scalar, 
hence most of its eigenstates will be a mixture of irreps of G, yet relation (jSJ) ensures that a subset 
of its eigenstates |[/ijv](E )A), are solvable and have good quantum numbers under the chain (jSJ). An 
Hamiltonian with partial dynamical symmetry is obtained by adding to H' the dynamical symmetry 
Hamiltonian, Hps 0, still preserving the solvability of states with (E) = (E ) 



H 



PDS 



Hns + H' 



(10) 



If the operators Tp = T^k^x span the entire irrep (a) of G, then the annihilation condition (|Sj) is 
satisfied for all A-states in (So), if none of the G irreps (S) contained in the Gdyn irrep [h N _ n ] belongs 
to the G Kronecker product (a) x (So). So the problem of finding interactions that preserve solvability 
for part of the states ([6]) is reduced to carrying out a Kronecker product. In this case, although the 
generators gi of G do not commute with H', their commutator does vanish when it acts on the solvable 
states (ED 



9i,H'\ ± 0, 
9i ,H' |[M(So)A) = 0, 9l eG. 



Eq. flllbj) follows from 



9i 



fit. 



ft 



9i,Tp 



+ 



Tp and the fact that 



(11a) 
(lib) 

involves 



a linear combination of G-tensor operators which satisfy Eq. (jSJ). The arguments for choosing the 
special irrep (S) = (S ) in Eq. (jSJ), which contains the solvable states, are based on physical grounds. 
A frequently encountered choice is the irrep which contains the ground state of the system. 

The above algorithm for constructing Hamiltonians with PDS of type I is applicable to any semisim- 
ple group. It can also address more general scenarios, in which relation (jSj) holds only for some states A 
in the irrep (S ) and/or some components A of the tensor T^ hn ^\. In this case, the Kronecker product 
rule does not apply, yet the PDS Hamiltonian is still of the form as in Eqs. (!9|)- (|T0|) . but now the solvable 
states span only part of the corresponding G-irrep. This is not the case in the quasi-exactly solvable 
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Hamiltonians, introduced in [23], where the solvable states form complete representations. The coexis- 
tence of solvable and unsolvable states, together with the availability of an algorithm, distinguish the 
notion of PDS from the notion of accidental degeneracy [26], where all levels are arranged in degenerate 
multiplets. 

An Hamiltonian with PDS of type I does not have good symmetry but some of its eigenstates do. 
The symmetry of the latter states does not follow from invariance properties of the Hamiltonian. This 
situation is opposite to that encountered in a spontaneous symmetry breaking, where the Hamiltonian 
respects the symmetry but its ground state breaks it. The notion of PDS differs also from the notion 
of quasi-dynamical symmetry [27J. The latter corresponds to a situation in which selected states in 
a system continue to exhibit characteristic properties (e.g., energy and B(E2) ratios) of a dynamical 
symmetry, in the face of strong symmetry-breaking interactions. Such an "apparent" persistence of 
symmetry is due to the coherent nature of the mixing in the wave functions of these states. In contrast, 
in a PDS of type I, although the symmetry is broken (even strongly) in most states, the subset of 
solvable states preserve it exactly. In this sense, the symmetry is partial but exact!. 

In what follows we present concrete constructions of Hamiltonians with PDS associated with the 
three DS chains (j2J) of the IBM. The partial symmetries in question involve continuous Lie groups. PDS 
can, however, be associated also with discrete groups which are relevant, e.g., to molecular physics. An 
example of a partial symmetry which involves point groups was presented in [28], employing a model of 
coupled anharmonic oscillators to describe the molecule XYq. The partial symmetry of the Hamiltonian 
allowed a derivation of analytic expressions for the energies and eigenstates of a set of unique levels. 
Furthermore, the numerical calculations required to obtain the energies of the remaining (non-unique) 
levels were greatly simplified since the Hamiltonian could be diagonalized in a much smaller space. 

2.1 U(5) PDS (type I) 

The U(5) DS chain of the IBM and related quantum numbers are given by [14J 

U(6) D U(5) D 0(5) D 0(3) 
ill | , (12) 

[N] (n d ) (t) n A L 

where the generators of the above groups are defined in Table dS] of the Appendix. For a given U(6) 
irrep [N], the allowed U(5) and 0(5) irreps are n d — 0, 1, 2, . . . , N and r = n d , rid — 2, ... or 1, respec- 
tively. The values of L contained in the 0(5) irrep (r), are obtained by partitioning r = 3n& + A, with 
tia, A > integers, and L = 2A, 2A — 2, 2A — 3, . . . , A. The multiplicity label n& in the 0(5) D 0(3) 
reduction, counts the maximum number of d-boson triplets coupled to L = [2S]. The eigenstates 
\[N](nd)(r)n&LM) are obtained with a Hamiltonian with U(5) DS which, for one- and two-body inter- 
actions, can be transcribed in the form 

#ds = en d + Ah d (n d + 4:) + BCo {5 ) + CCo(3) • (13) 

Here fid and n d (n d + 4) are the linear and quadratic Casimir operators of U(5), respectively, and Cq 
denotes the quadratic Casimir operator of G, as defined in the Appendix. The Casimir operators of U(6) 
are omitted from Eq. ( fl3|) since they are functions of the total boson number operator, N = n s + n d , 
which is a constant for all iV-boson states. The spectrum of H-ds is completely solvable with eigenenergies 

£ds = 6n d + An d (n d + 4) + BT(T + 3)+ CL(L + 1). (14) 

The U(5)-DS spectrum of Eq. (fl4"|) resembles that of an anharmonic spherical vibrator, describing 
quadrupole excitations of a spherical shape. The splitting of states in a given U(5) multiplet, (n d ), 
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Table 1: Normalized one- and two-boson U(5) tensors. 
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is governed by the 0(5) and 0(3) terms in ifrjs f[T3"|) . The lowest U(5) multiplets involve states with 
quantum numbers (n d = 0, r = 0, L = 0), (n^ = 1, r = 1, L = 2), and (n d = 2, r = 0, L = 0), 
(n d = 2, r = 2, L = 2,4). 

The construction of Hamiltonians with U(5)-PDS is based on identification of n-boson operators 
which annihilate all states in a given U(5) irrep (n d ). A physically relevant choice is the irrep rid — 
which consists of the ground state, with r = L = 0, built of iV s-bosons 



\[N],n d 



T 



L = 0) = (N\)- 1/2 (sY\ty ■ 



(15) 



Considering U(5) tensors, B^ n ^ n ^^ n ^ lm , composed of n bosons of which rid are <i-bosons then, clearly, 



the Hermitian conjugate of such operators with rid 7^ will annihilate the n d = state of Eq. (115]) . 
Explicit expressions for n-boson U(5) tensors, with n = 1,2 are shown in Table [TJ From them one can 
construct the following one- and two-body Hamiltonian with U(5)-PDS 



Hpds = e d d ] ■ d + u 2 s ] d ] ■ ds + v 2 s ] S ■ {dd) {2) + H .c. 
+ c L(d j d^ ■ (dd)W , 

L=0,2,4 

where H.c. means Hermitian conjugate. By construction, 

H PDS \[N},n d = r = L = 0) = 0. 
Using Eq. (11821) . we can rewrite Hp us 111 the form 

Hpos — H DS + V 2 , 
where Hp>s is the U(5) dynamical symmetry Hamiltonian, Eq. fTl3|) . and V 2 is given by 



(16) 



(17) 



V 2 



v 2 



sW ■ (dd)W + h.c. 1 = v 2 n( 2 ) • = v 2 ■ 11® . 



(19) 



The operators and Up are defined in Eq. ([TJ). The V 2 term breaks the U(5) DS, however, it still 
has the U(5) basis states with (n d — r — L — 0), Eq. (EETj), and (rid = T = L = 3) as zero-energy 
eigenstates 



V 2 \[N],n d 



T 



3) 







(20) 
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The last property follows from the U(5) selection rules of V2, An d = ±1, and the fact that the irreps 
(n d = 4, r = 0, 2, 4) and (rid — 2, r = 0, 2) do not contain an L = 3 state. Altogether, Hpos ( fl8|) is not 
diagonal in the U(5) chain, but retains the following solvable U(5) basis states with known eigenvalues 

\[N],n d = r = L = 0) E PDS = 0, (21a) 

\[N},n d = r = L = 3) E PDS = 3e + 21 A + 18B + 12C . (21b) 

As will be discussed in Section El this class of Hamiltonians with U(5)-PDS of type I is relevant to 
first-order quantum shape-phase transitions in nuclei. 

A second class of Hamiltonians with U(5)-PDS can be obtained by considering the interaction 



V = v 



(s ] ) 2 d-d + H.c. 



(22) 



This interaction breaks the U(5) DS, however, it still has selected U(5) basis states as zero-energy 
eigenstates 

V \[N},n d = r = N,L) = 0, (23a) 

V \[N],n d = T = N-l,L) = 0, (23b) 

where L takes the values compatible with the 0(5) D 0(3) reduction. These properties follow from 
the fact that s 2 annihilates states with n s — 0, 1 (n d = N, N — 1) and d ■ d annihilates states with 
n d = t [14J. Adding the interaction Vq to the U(5) dynamical symmetry Hamiltonian Hps (113j) . we 
obtain the following Hamiltonian with U(5)-PDS 

H'pds = H DS + V . (24) 

H'pds i s n °t diagonal in the U(5) chain, but retains the following solvable U(5) basis states with known 
eigenvalues 

\[N},n d = r = N,L) E' PDS = N[e d + A(N + 4) + B(N + 3)] + CL(L + l) , (25a) 

\[N},n d = r = N -1,L) E' PDS = (N-l)[e d + A(N + 3)+B(N + 2)}+CL(L + l) . (25b) 

The Hamiltonian H' PDS ( 124"]) with U(5)-PDS of type I, contains terms from both the U(5) and 0(6) 
chains (j2J), hence preserves the common segment of subalgebras, 0(5) D 0(3). As such, it exhibits 
also an 0(5)-PDS of type II, to be discussed in Sectional As will be shown in Section |H1 this class of 
Hamiltonians is relevant to second-order quantum shape-phase transitions in nuclei. 

2.2 SU(3) PDS (type I) 

The SU(3) DS chain of the IBM and related quantum numbers are given by [15] 



U(6) D SU(3) D 0(3) 
I I I , (26) 

[N] (A,/i) K L 

where the generators of the above groups G are defined in Table [TBI of the Appendix. For a given U(6) 
irrep [N], the allowed SU(3) irreps are (A,/i) = (2N — 4k — 6m, 2k) with k,m non-negative integers, 
such that, A,/i > 0. The multiplicity label K is needed for complete classification and corresponds geo- 
metrically to the projection of the angular momentum on the symmetry axis. The values of L contained 
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in the above SU(3) irreps are L — K, K + 1, K + 2, . . . , K + max{A, //}, where K — 0, 2, ... , min{A, fi}; 
with the exception of K = for which L = 0, 2, . . . , max{A, /i}. The states |[iV](A, [i)KLM) form 
the (non-orthogonal) Elliott basis |4J and the Vergados basis |[iV](A, [i)xLM ) [30J is obtained from it 
by a standard orthogonalization procedure. The two bases coincide in the large-N limit and both are 
eigenstates of a Hamiltonian with SU(3) DS. The latter, for one- and two-body interactions, can be 
transcribed in the form 



Hds = h<2 



4u ( 3)+2iV(2iV + 3) +CC oi 3) , (27) 



where Cq is the quadratic Casimir operator of G, as defined in the Appendix. The spectrum of H DS is 
completely solvable with eigenenergies 

£ds = h 2 [-f 2 (\,fi) + 2N(2N + 3)}+CL(L + l) 

= h 2 6 [2N(k + 2m) - k(2k - 1) - 3m(2m - 1) - 6km] + CL(L + 1) , (28) 

where f 2 (X, /i) = A 2 + fi 2 + Xfi + 3A + 3/i and (A, fi) = {2N — 4k — 6m, 2k). The spectrum resembles that 
of an axially-deformed rotovibrator and the corresponding eigenstates are arranged in SU(3) multiplets. 
In a given SU(3) irrep (A,/x), each K- value is associated with a rotational band and states with the 
same L, in different K-bands, are degenerate. The lowest SU(3) irrep is (2N, 0), which describes the 
ground band g(K = 0) of a prolate deformed nucleus. The first excited SU(3) irrep (2N — 4, 2) contains 
both the (3(K = 0) and j(K = 2) bands. States in these bands with the same angular momentum 
are degenerate. This /3-7 degeneracy is a characteristic feature of the SU(3) limit of the IBM which, 
however, is not commonly observed [31]. In most deformed nuclei the /3 band lies above the 7 band. In 
the IBM framework, with at most two-body interactions, one is therefore compelled to break SU(3) in 
order to conform with the experimental data. To do so, the usual approach has been to include in the 
Hamiltonian terms from other chains so as to lift the undesired ^-7 degeneracy. Such an approach was 
taken by Warner Casten and Davidson (WCD) in [32], where an 0(6) term was added to the SU(3) 
Hamiltonian. However, in this procedure, the SU(3) symmetry is completely broken, all eigenstates 
are mixed and no analytic solutions are retained. Similar statements apply to the description in the 
consistent Q formalism (CQF) [33] , where the Hamiltonian involves a non-SU(3) quadrupole operator. 
In contrast, partial SU(3) symmetry, to be discussed below, corresponds to breaking SU(3), but in a 
very particular way so that part of the states (but not all) will still be solvable with good symmetry. As 
such, the virtues of a dynamical symmetry (e.g., solvability) are fulfilled but by only a subset of states. 

The construction of Hamiltonians with SU(3)-PDS is based on identification of n-boson operators 
which annihilate all states in a given SU(3) irrep (A, //), chosen here to be the ground band irrep (2N, 0). 
For that purpose, we consider the following boson-pair operators with angular momentum L = 0, 2 

pt = d ] . d ] _ 2 ( s t)2 ) (29a) 
Pi = 2S^ + V7(dU^. (29b) 

As seen from Table [2, these operators are proportional to two-boson SU(3) tensors, -SL(A M )^ m > with 
n = 2 and (A,/i) = (0,2) 

^[2](o,2)o;oo = '^/f^o ' (30a) 
B U (0,2)0^ = ^/f^V ( 30b ) 
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Table 2: Normalized one- and two-boson SU(3) tensors. For the indicated irreps the labels of the Vergados 
basis (x) and Elliott basis (K) are identical, \ = K. 
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The corresponding Hermitian conjugate boson-pair annihilation operators, Pq and Pip,-, transform as 
(2,0) under SU(3), and satisfy 

P \[N](2N,0)K = 0,LM) = 0, 

P 2li \[N](2N,0)K = 0,LM) = 0. (31) 
Equivalently, these operators satisfy 

P |c; AO = 

P 2 p\c;N) = (32) 

where 

\c;N) = (N\)-V 2 (bi) N \0) , 6+ = (V2 dj + st)/V3 . (33) 

The state |c; AT) is a condensate of A" bosons and is obtained by substituting the SU(3) equilibrium 
deformations in the coherent state of Eq. (j3J), |c; iV) = |/3 = a/2, 7 = 0; A 7 }. It is the lowest-weight state 
in the SU(3) irrep (A, /x) = (2N, 0) and serves as an intrinsic state for the SU(3) ground band [33]. The 
rotational members of the band |[A / ](2A r , 0)K = 0, LM), Eq. fj3T|) . are obtained by angular momentum 
projection from |c; AT). The relations in Eqs. f[3~Tj) - f[3"2"j) follow from the fact that the action of the 
operators leads to a state with N — 2 bosons in the U(6) irrep [N — 2], which does not contain the 
SU(3) irreps obtained from the product (2, 0) x (2N, 0) = (2N + 2, 0) © (2N, 1) © (2N - 2, 2). 

Following the general algorithm, a two-body Hamiltonian with partial SU(3) symmetry can now be 
constructed as [231135] 

H(h ,h 2 ) = h P^P Q + h 2 P^P 2 , (34) 

where P 2 ^ = (— ) fl P 2 ~p- For h 2 = h , the Hamiltonian is an SU(3) scalar, related to the quadratic 
Casimir operator of SU(3) 

H(h = h 2 ) = h 2 [pJPo + pt.pJ =h 2 L(7 SU( 3 )+ 2A>(2A> + 3)1 . (35) 
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For h = —5h 2 , the Hamiltonian transforms as a (2,2) SU(3) tensor component. For arbitrary h ,h 2 
coefficients, H(ho,h 2 ) is not an SU(3) scalar, nevertheless, the relations in Eqs. (13T1) - (132"1) ensure that 
it has a solvable zero-energy band with good SU(3) quantum numbers (A,/i) = (27V, 0). When the 
coefficients h , h 2 are positive, H(h , h 2 ) becomes positive definite by construction, and the solvable 
states form its SU(3) ground band. 

H(h , h 2 ) of Eq. (j3"4"|) has additional solvable eigenstates with good SU(3) character. This comes 
about from Eq. (1321) and the following properties 

P L Jc;N) = 



= 5 Lt2 5^ 2 24P{ 2 , L = 0,2 . (36) 

These relations imply that the sequence of states 

\k) oc (p{ 2 ) k \ C] N-2k) , (37) 

are eigenstates of H(ho, h 2 ) with eigenvalues Ek = 6h 2 (2N + 1 — 2k) k. A comparison with Eq. (1281) 
shows that these energies are the SU(3) eigenvalues of H(h = h 2 ), Eq. (1331) . and identify the states \k) 
to be in the SU(3) irreps (2N — 4k, 2k). They can be further shown to be the lowest-weight states in 
these irreps. Furthermore, Pq satisfies 

Po\k) = 0, (38) 

or equivalently, 

P Q \[N]{2N -4k,2k)K = 2k,LM) = 0. (39) 

The states \k) (1371) are deformed and serve as intrinsic states representing 7 fc bands with angular 
momentum projection (K = 2k) along the symmetry axis [31]. In particular, as noted earlier, \k = 
0) = \c; N) represents the ground band (K = 0) and \k = 1) is the 7-band (K = 2). The rotational 
members of these bands, \[N](2N - 4k, 2k)K = 2k, LM) oc V LM \k), Eq. (Pj), can be obtained by 0(3) 
projection from the corresponding intrinsic states \k). Relations (1351) and (13il are equivalent, since the 
angular momentum projection operator, Vlm, is composed of 0(3) generators, hence commutes with 
Pq. The intrinsic states \k) break the 0(3) symmetry, but since the Hamiltonian in Eq. (1341) is 0(3) 
invariant, the projected states with good L are also eigenstates of H(ho, h 2 ) with energy Ek and with 
good SU(3) symmetry. For the ground band {k = 0) the projected states span the entire SU(3) irrep 
(2N, 0). For excited bands {k 7^ 0), the projected states span only part of the corresponding SU(3) 
irreps. There are other states originally in these irreps (as well as in other irreps) which do not preserve 
the SU(3) symmetry and therefore get mixed. In particular, the ground g(K = 0), and ^(K = 2) bands 
retain their SU(3) character (2N, 0) and (27V — 4, 2) respectively, but the first excited K — band 
is mixed. This situation corresponds precisely to that of partial SU(3) symmetry. An Hamiltonian 
H(h , h 2 ) which is not an SU(3) scalar has a subset of solvable eigenstates which continue to have good 
SU(3) symmetry. 

All of the above discussion is applicable also to the case when we add to the Hamiltonian of Eq. (1341) 
the Casimir operator of 0(3), and by doing so, convert the partial SU(3) symmetry into partial dynam- 
ical SU(3) symmetry. The additional rotational term contributes just an L(L + 1) splitting but does 
not affect the wave functions. The most general one- and two-body Hamiltonian with SU(3)-PDS can 
thus be written in the form 

H PDS = H(h ,h 2 ) + CC o(3) = H DS + (h -h 2 )P^P . (40) 



L,fi j 



P t 
r 2,2 



\c;N) = 5 L , 2 5 Mi2 6(27V + 3)|c;7V) 



Pl,h 1 P' 



2,2 



P, 



2.2 
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Here Hds is the SU(3) dynamical symmetry Hamiltonian, Eq. (12 7p. with parameters h% and C. The 
PqPq term in Eq. ( l4"0j) is not diagonal in the SU(3) chain ( )26|) . but Eqs. ( )3T|) and (|39|) ensure that it 
annihilates a subset of states with good SU(3) quantum numbers. Consequently, Hpds retains selected 
solvable bands with good SU(3) symmetry. Specifically, the solvable states are members of the ground 
g(K = 0) band 

\N, (2N,0)K = 0,L) L = 0,2,4, ... , 2N (41a) 
E PDS = CL{L + 1) (41b) 

and 7 fc (K = 2k) bands 

\N,(2N -Ak,2k)K = 2k,L) L = K, K + 1, . . . , (2N - 2k) (42a) 
E PDS = h 2 6k(2N -2k + l) + CL(L + 1) k > . (42b) 

The solvable states f l4Tj) -f T42|) are those projected from the intrinsic states \k) of Eq. fl37|) . and are 
simply selected members of the Elliott basis <Pe{{\ fi)KLM) [I]. In particular, the states belonging 
to the ground and 7 bands are the Elliott states <M(2iV, 0)K = 0,LM) and 4> E ((2N - A,2)K = 
2, LM) respectively. Their wave functions can be expressed in terms of the orthonormal Vergados 
basis, \Ev((A, n)xLM) [30J. For the ground band and for members of the 7 band with L odd, the 
Vergados and Elliott bases are identical. The Elliott states in the 7 (if = 2) band with L even are 
mixtures of Vergados states in the SU(3) irrep (2N — 4, 2) 



</>e((2iV — 4, 2)K = 2, LM) 



l m v ((2N - 4, 2)x = 2, LM) - arg? ^ V ((2N - 4, 2)x = 0, LM) j /x£> , (43) 

where x^, are known coefficients which appear in the transformation between the two bases [30] . 

Since the wave functions of the solvable states are known, it is possible to obtain analytic expressions 
for matrix elements of observables between them. For calculating E2 rates, it is convenient to rewrite 
the relevant E2 operator, Eq. ([3]), in the form 

T{E2) = aQ {2) + 6Ii {2) , (44) 

where Q (2) is the quadrupole SU(3) generator [Q^ = - (y/7/2) U^] and n^ 2 ) is a (2,2) tensor 
under SU(3). The B(E2) values for g — > g and g — > 7 transitions 



B(E2;g,L-*g,L') = 

I (<M(2iV, 0)K = 0, V) \\aQW + 9 IT 2 ) \\(f> E ((2N, 0)K = 0, L))\ 2 

(2L + 1) 

B{E2; 1 ,L^g,L') = 

]2 I (Mm 0)K = 0, V) I |H( 2 ) \\M{ZN - 4, 2)K = 2, L))\ 2 



(45a) 



e ^ifff — • (45b) 

can be expressed in terms of E2 matrix elements in the Vergados basis, for which analytic expressions 
are available [T5|l36]. 

The Hamiltonian of Eq. (f40|) . with SU(3)-PDS, was used in [35] to describe spectroscopic data of 
168 Er. The experimental spectra [32] of the ground g(K = 0^), ^{K = 2f) and K = 0^ bands in this 
nucleus is shown in Fig. [21 and compared with an exact DS, PDS and broken SU(3) calculations. 
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Figure 2: Spectra of 168 Er (N = 16). Experimental energies (EXP) are compared with IBM calculations in an 
exact SU(3) dynamical symmetry [SU(3)], in a broken SU(3) symmetry (WCD) |32j and in a partial dynamical 
SU(3) symmetry (PDS). The latter employs the Hamiltonian of Eq. (|4U|) with /iq = 8, h 2 = 4, C = 13 keV. 
Adapted from [35] . 



According to the previous discussion, the SU(3)-PDS spectrum of the solvable ground and 7 bands is 
given by 

E g (L) = CL{L + 1) , 

E 7 (L) = 6h 2 (2N - 1) + CL(L + 1) . (46) 

Hpds P0|) is specified by three parameters (N=16 for 168 Er according to the usual boson counting). 
The values of C and h 2 were extracted from the experimental energy differences [E(2+) — E(0^)} and 
[E{2^) — E(2+)] respectively. The spectrum of an exact SU(3) DS, Eq. fT25|) . obtained for h = h 2 , 
deviates considerably from the experimental data since empirically the K = and j(K = 2^) bands 
are not degenerate. In the SU(3)-PDS calculation, the parameter h Q was varied so as to reproduce 
the bandhead energy of the K = band. The prediction for other rotational members of the K = 
0^, 0^", l\ bands is shown in Fig. |2j Clearly, the SU(3) PDS spectrum is an improvement over the 
schematic, exact SU(3) dynamical symmetry description, since the /3-7 degeneracy is lifted. The ground 
and 7 bands are solvable, still, the quality of the calculated PDS spectrum is similar to that obtained 
in the broken-SU(3) calculation (WCD) [32]. 

The resulting SU(3) decomposition of the lowest bands for the SU(3)-PDS calculation is shown in 
Fig. |3j and compared [37] to the conventional broken-SU(3) calculations WCD [32] and CQF [33J. In 
the latter calculations all states are mixed with respect to SU(3). In contrast, in the PDS calculation, 
the good SU(3) character, (32, 0) for the ground band and (28, 2) for 7 band, is retained, while the 
K = O2 band contains about 13% admixtures into the dominant (28, 2) irrep. Thus, the breaking of 
the SU(3) symmetry induced by Hpos P0|) is not small, and can lead to an appreciable SU(3) mixing 
in the remaining non-solvable states. Nevertheless, the special solvable states carry good SU(3) labels. 
The SU(3) symmetry is therefore partial but exact. 

Electromagnetic transitions are a more sensitive probe to the structure of states, hence are an 



15 



£ 80 

I 60 

i 40 
o 

fa 20 



■ PDS 


CQF 


_ wcd : 




■ 





100 



S 60 



5 40 



ft 20 



100 



!J 40 

- 



I 20 



(32,0) (24,4) (32,0) (24,4) (32,0) (24,4) 
(28,2) (28,2) (28,2) 




(28,2) (20,6) (28,2) (20,6) (28,2) (20,6) 
(24,4) (24,4) (24,4) 





PDS 


CQF 


WCD 












K=0, 






-■ 


ill .. 




L 



(32,0) (26,0) (32,0) (26,0) (32,0) (26,0) 
(28,2) (20,6) (28,2) (20,6) (28,2) (20,6) 
(24,4) (16,8) (24,4) (16,8) (24,4) (16,8) 



Figure 3: SU(3) decomposition of wave functions of the ground (K = Oi), 7 (K = 2i), and K = O2 bands 
of 168 Er (N = 16) in a SU(3)-PDS calculation and in broken-SU(3) calculations: CQF [33J and WCD [32]. 
Adapted from [37] . 



important indicator for testing the predictions of SU(3)-PDS. Based on the analytic expressions of 
Eq. (Ho]) , the parameters a and 9 of the E2 operator can be extracted from the experimental values of 
B(E2; 0+ ->• 2+) and B(E2; 0+ ->• 2+). The corresponding ratio for 168 Er is 9/a = 4.261. As shown in 
Table |3l the resulting SU(3) PDS E2 rates for transitions originating within the 7 band are found to 
be in excellent agreement with experiment and are similar to the WCD broken-SU(3) calculation [32J. 
In particular, the SU(3) PDS calculation reproduces correctly the (7 — > —> g) strengths. The 

only significant discrepancy is that for the 8^ — > 7^ transition which is very weak experimentally, with 
an intensity error of 50% and an unknown Ml component [32] ■ The expressions in Eq. (145]) imply 
that 7 — > g B(E2) ratios are independent of both E2 parameters a and 9. Furthermore, since the 
ground and 7 bands have pure SU(3) character, (2N, 0) and (2A^ — 4, 2) respectively, the corresponding 
wave-functions do not depend on parameters of the Hamiltonian and hence are determined solely by 
symmetry. Consequently, the B(E2) ratios for 7 — > g transitions quoted in Table [3] are parameter-free 
predictions of SU(3) PDS. The agreement between these predictions and the data confirms the relevance 
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Table 3: B(E2) branching ratios from states in the 7 band in 168 Er. The column EXP lists the experimental 
ratios [S2], PDS is the SU(3) partial dynamical symmetry calculation [35J and WCD is the broken SU(3) 
calculation [32]. Adapted from [35] , 
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of SU(3)-PDS to the spectroscopy of 168 Er. 

PDS has important implications not only for the structure of the pure (solvable) states, but it 
also affects the mixing of the remaining (non-solvable) states. Of particular interest is the nature of 
the lowest K = 2 excitation, for which the role of double-7-phonon admixtures is still subject to 
controversy [38J. A closer look at the SU(3) decomposition in Fig. [31 shows that in the SU(3)-PDS 
calculation, the K = O2 band is mixed and has the structure 



\L,K = 2 ) = A 1 M(2N-4,2)K = 0,L) + A 2 4> E ((2N-8,4:)K = 0,L) 

+A 3 4> E ((2N-6,0)K = 0,L) . (47) 

Here (pE denote states orthogonal to the solvable 7^ =2 fc Elliott's states. The notation K = signifies 
that K = is only the dominant component in these states. For example, 

M(2N-A,2)K = 0,L) = 

W((2N - 4, 2)x = 0, L) + 4oW((2iV - 4, 2)x = 2, L)J /x%> , (48) 

is the state orthogonal to the Elliott state in Eq. f )43|) . For 168 Er, the K = O2 band was found [37] to 
contain 9.6% (26,0) and 2.9% (24,4) admixtures into the dominant (28,2) irrep. Using the geometric 
analogs of the SU(3) bands [39J, (2N - A,2)K = ~ /3, (2N - 8,4)K = ~ (V2/3 2 + 7 2 x=o), 
(2N — 6, 0)K = ~ (/3 2 — ^7^=0); the wave function of Eq. fl4T|) can be expressed in terms of the 
probability amplitudes for single- and double- phonon K = excitations 

A/3 = A 1 , A 7 2 = (A 2 - V2A 3 )/V3 , A 02 = (V2A 2 + A 3 )/Vs . (49) 
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Table 4: Comparison of calculated (Calc) and experimental (Exp) [12] absolute B(E2) values [W.u.] for 
transitions from the K = O2 band in 168 Er. PDS is the SU(3) partial dynamical symmetry calculation [37], 
while WCD [32] and CQF [33J are broken SU(3) calculations. Adapted from [37]. 



Exp 

Transition B(E2) Range 


Calc 

PDS WCD CQF 


2£ =02 -> 0+ 0.4 0.06-0.94 
2£ =02 ->• 2+ 0.5 0.07-1.27 
2+ =0a ^4+ 2.2 0.4-5.1 
2i =02 ->• 2+ a) 6.2 (3.1) 1-15 (0.5-7.5) 
2+ =02 ->• 3+ fl ) 7.2 (3.6) 1-19 (0.5-9.5) 


0.65 0.15 0.03 
1.02 0.24 0.03 
2.27 0.50 0.10 
4.08 4.2 4.53 
7.52 7.9 12.64 



a ) The two numbers in each entry correspond to an assumption of pure E2 
and (in parenthesis) 50% E2 multipolarity. 



It follows that in the PDS calculation, the K = O2 band in 168 Er contains admixtures of 12.4% 7^ =0 
and 0.1% (3 2 into the /3 mode, i.e. 12.5% double-phonon admixtures into the dominant single-phonon 
component. These findings support the conventional single-phonon interpretation for this band with 
small but significant double-7-phonon admixture. 

General properties of the K = 2 band have been studied [37] by examining the SU(3)-PDS 
Hamiltonian of Eq. (134"1) . The empirical value of the ratio of K = O2 and 7 bandhead energies 
E(0 2 )/[E(2+) — E(2g)] = 0.8 — 1.8, in the rare-earth region [4"0~||4T] constrains the parameters of 

H(ho, h 2 ) to be in the range 0.7 < h /h 2 < 2.4. In general, one finds that the K = 2 wave function 
retains the form as in Eq. f|47|) and, therefore, a three-band mixing calculation is sufficient to describe 
its structure. The SU(3) breaking and double-phonon admixture is more pronounced when the K = 2 
band is above the 7 band. For most of the relevant range of ho/h 2 , corresponding to bandhead ratio 
in the range 0.8 — 1.65, the double-phonon admixture is at most ~ 15%. Only for higher values of the 
bandhead ratio can one obtain larger admixtures and even dominance of the 7x =0 component in the 
K = O2 wave function. 

An important clue to the structure of K = 2 collective excitations comes from E2 transitions 
connecting the K = 2 and K = 0i bands. If we recall that only the ground band has the SU(3) 
component (A,//) = (2N,0), that Q^ 2 \ as a generator, cannot connect different SU(3) irreps and that 
(2,2) tensor under SU(3), can connect the (2N,0) irrep only with the (2N — 4,2) irrep, we 
obtain the following expression for the B(E2) values of K = 2 — > g transitions 



B(E2;K = 2 ,L^ g,L') 

a2 q2 1 (<M(2iv, o)k = 0, V) I |n( 2 ) \\M(2N - 4, 2)K = 0, L))\ z 
A e 6 (2T+T) • ! 



Employing Eq. ( 1481) . this B(E2) value can be expressed in terms of known B(E2) values in the Vergados 
basis [15j[36]. Using Eq. (145b|) . the E2 parameter 9 can be determined from the known 2+ — > 0+ E2 
rates, and for 168 Er is found to be 9 2 = 2.175 W.u. As seen from Eq. (150]) . the B(E2) values for 
K = O2 — > g transitions are proportional to {Ap) 2 , hence, they provide a direct way for extracting the 
amount of SU(3) breaking and the admixture of double-phonon excitations in the K = 2 wave function. 
In Table H] we compare the predictions of the PDS and broken-SU(3) calculations with the B(E2) values 
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deduced from a lifetime measurement of the 2^ =0a level in 168 Er [42j (^ ne indicated range for the B(E2) 
values correspond to different assumptions on the feeding of the level). The PDS and WCD calculations 
are seen to agree well with the empirical values but the CQF calculation under predicts the measured 
K = 2 — > g data. 

The SU(3)-PDS discussed above, is relevant to rotational-vibrational states of a prolate deformed 
shape, with equilibrium deformations (f3 = a/2, 7 = 0) and a symmetry z-axis. It is also possible to 
identify an SU(3)-PDS corresponding to an oblate shape with equilibrium deformations ((3 = a/2,7 = 
7r/3) and a symmetry y-axis. The Hamiltonian with SU(3)-PDS has the same form as in Eq. fj4*0]) but 
the L = 0, 2 boson-pairs are obtained from Eq. (1291) by a change of phase: — > —s\ s — > —s. The 
generators and quadratic Casimir operator of SU(3) are listed in the Appendix. The relevant SU(3) 
irreps, (A, p) = (2k, 2N — 4k — 6m), are conjugate to the SU(3) irreps, (A, /i), encountered in the SU(3) 
chain, Eq. ([26]). 

2.3 0(6) PDS (type I) 

The 0(6) DS chain of the IBM and related quantum numbers are given by [TB] 

U(6) D 0(6) D 0(5) D 0(3) 
ill | , (51) 

[N] (a) (r) n A L 

where the generators of the above groups are listed in Table [TBI of the Appendix. For a given U(6) irrep 
[N], the allowed 0(6) and 0(5) irreps are a = N, N — 2, ... or 1, and r = 0, 1, ... a, respectively. The 
0(5) D 0(3) reduction is the same as in the U(5) chain. The eigenstates \ [N](a) (r)n A LM) are obtained 
with a Hamiltonian with 0(6) DS which, for one- and two-body interactions, can be transcribed in the 
form 



-C 0( 6) + N(N + 4) + B C 0(5 ) + C C 0( 3) • (52) 



The quadratic Casimir operators, Cg, are defined in the Appendix. The spectrum of _£/ds is completely 
solvable with eigenenergies 

£ds = h (N-(r)(N + (T + 4) + BT(T + 3)+ C L(L + 1) . (53) 

The spectrum resembles that of a 7- unstable deformed rotovibrator, where states are arranged in 0(6) 
multiplets with quantum number a. The ground band corresponds to the 0(6) irrep with a = N. 
The splitting of states in a given 0(6) multiplet is governed by the 0(5) and 0(3) terms in H DS (152"|) . 
The lowest members in each band have quantum numbers (r = 0, L = 0), (r = 1, L = 2) and 
(r = 2, L = 2,4). 

The construction of Hamiltonians with 0(6)-PDS is based on identification of n-boson operators 
which annihilate all states in a given 0(6) irrep, (a), chosen here to be the ground band irrep (er) = (N). 
For that purpose, a relevant operator to consider is 



pt = £ ■$ - (s f ) 2 . (54) 

t 



As seen from Table [SJ the above operator is proportional to a two-boson 0(6) tensor, B^ n ^ a ^ r 
with n = 2 and o = r = L = 



-^[2]<0}(0)0;00 — 2V3^ ' 
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Table 5: Normalized one- and two-boson 0(6) tensors. 



n 


cr 


r 




f 


-i 
l 


1 











1 


1 


1 





2 


9 
Z 


o 
z 


u 


n 
u 


n 
u 


2 


2 


1 





2 


2 


2 


2 





2 


2 


2 


2 





4 


2 















5, 



t 

n] (a) (t)tia ;&ra 



S 1 

rn 



(-1) 



5 {<#S)f 



The corresponding Hermitian conjugate boson-pair annihilation operator, Pq, transforms also as (a) = 
(0) under 0(6) and satisfies 

P \[N)(N)(r)n A LM) =0 . (56) 

Equivalently, this operator satisfies 

P |c; AO = (57) 



where 



\c;N) = (Niy^ibl^lO) 

b\ = [cos7 4 + sin7(4 + ^ t _ 2 )/v / 2 + s t ]/v / 2 . (5? 



The state \c; N) is obtained by substituting the 0(6) equilibrium deformations in the coherent state of 
Eq. ([3]), |c; N) = \(3 = l,j;N). It is the lowest-weight state in the 0(6) irrep (cr) = (N) and serves as an 
intrinsic state for the 0(6) ground band [13] • The rotational members of the band, \[N]{N)(T)riALM) , 
Eq. (156]) . are obtained by 0(5) projection from |c; N) and span the entire 0(6) irrep (cr) = (N). The 
relations in Eqs. (l56l) - fl5"Tj) follow from the fact that the action of the operator P leads to a state with 
N — 2 bosons in the U(6) irrep [N — 2], which does not contain the 0(6) irrep (N) obtained from the 
product of (0) x (N). 

Since both Pq and Pq (|5^|) are 0(6) scalars, they give rise to the following 0(6)-invariant interaction 



C oie) +N(N + A) , (59) 



which is simply the 0(6) term in Hds, Eq. (|52j) . with an exact 0(6) symmetry. Thus, in this case, 
unlike the situation encountered with SU(3)-PDS, the algorithm does not yield an 0(6)-PDS of type I 
with two-body interactions. In the IBM framework, an Hamiltonian with a genuine 0(6)-PDS of this 
class, requires higher-order terms. A construction of three-body Hamiltonians with such property will 
be presented in Subsection 18.21 
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3 PDS type II 



PDS of type II corresponds to a situation for which all the states of the system preserve part of the 
dynamical symmetry, 

G D Gi D G 2 D . . . D G n . (60) 

In this case, there are no analytic solutions, yet selected quantum numbers (of the conserved symmetries) 
are retained. This occurs, for example, when the Hamiltonian contains interaction terms from two 
different chains with a common symmetry subalgebra, e.g., 

G D | ^ | D G 2 D . . . D G n (61) 

If G\ and G' x are incompatible, i.e., do not commute, then their irreps are mixed in the eigenstates of 
the Hamiltonian. On the other hand, since G 2 and its subalgebras are common to both chains, then 
the labels of their irreps remain as good quantum numbers. An Hamiltonian based on a spectrum 
generating algebra Gq and a symmetry algebra G n common to all chains of Go, has by definition, a 
G n -PDS of type II, albeit a trivial one [e.g., G n = 0(3)]. Therefore, the notion of PDS type II is 
physically relevant when the common segment of the two DS chains, contains subalgebras which are 
different from the symmetry algebra, i.e., G 2 ^ G n in Eq. (IBTj) . 

An alternative situation where PDS of type II occurs is when the Hamiltonian preserves only some 
of the symmetries Gj in the chain (|60|) and only their irreps are unmixed. A systematic procedure for 
identifying interactions with such property was proposed in [44J. Let Gi D G 2 D G 3 be a set of nested 
algebras which may occur anywhere in the reduction ( I60p . in-between the spectrum generating algebra 
Go and the invariant symmetry algebra G n . The procedure is based on writing the Hamiltonian in terms 
of generators, gi, of G\, which do not belong to its subalgebra G 2 . By construction, such Hamiltonian 
preserves the G\ symmetry but, in general, not the G 2 symmetry, and hence will have the G\ labels as 
good quantum numbers but will mix different irreps of G 2 . The Hamiltonians can still conserve the G3 
labels e.g., by choosing it to be a scalar of G3. The procedure for constructing Hamiltonians with the 
above properties involves the identification of the tensor character under G 2 and G3 of the operators gi 
and their products, gigj . . . gu- The Hamiltonians obtained in this manner belong to the integrity basis 
of G3-scalar operators in the enveloping algebra of G\ and, hence, their existence is correlated with 
their order. In the discussion below we exemplify the two scenarios for constructing Hamiltonians with 
PDS of type II within the IBM framework. 

3.1 0(5) PDS (type II) 

An example of mixing two incompatible chains with a common symmetry subalgebra is the U(5) and 
0(6) chains in the IBM [IS] 

U(6) D { } D 0(5) D 0(3) . (62) 

The corresponding quantum numbers were discussed in Subsection 12.11 and 12.31 The most general 
Hamiltonian which conserves the 0(5) symmetry, involves a combination of terms from these two 
chains, and for one- and two-body interactions reads 

#0(5) = en d + a n d {n d + 4) + A P$P + B G 0(5 ) + C G 0(3 ) • (63) 
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The d-boson number operator, n^, is the linear Casimir of U(5), the 0(6)-pairing term, PqPq, is related 
to the Casimir operator of 0(6), Eq. (|59|) . and the Casimir operators of 0(5) and 0(3) can be replaced by 
their eigenvalues B r(r + 1) + C L(L+1). In this case, all eigenstates of Ho(B) have good 0(5) symmetry 
but, with a few exceptions, none of them have good U(5) symmetry nor good 0(6) symmetry, and hence 
only part of the dynamical symmetry of each chain in Eq. (1621) is observed. These are precisely the 
defining features of 0(5)-PDS of type II. In general, #0(5) dMD mixes U(5) irreps, characterized by n d , 
as well as mixes irreps of 0(6) characterized by a, but retains the 0(5) D 0(3) labels, (r, L) as good 
quantum numbers. The conserved 0(5) symmetry has important consequences which will be discussed 
briefly below [35] . 

The first four terms in ( 1631) are 0(5) scalars, hence level-spacings within an 0(5) irrep (r-multiplet) 
are the same throughout the U(5)-0(6) transition region. They are given only by CL(L + 1) and are 
independent of the values of the other parameters in Ho(5)- The various multiplets with the same 
value of t are distinguished by another label, v = 1, . . . which indicates its relative position in the 
spectrum. The actual positions of the various r-multiplets, as well as the wave functions of their states, 
are determined by diagonalization of the Hamiltonian (163 p . The wave function of any state in the z/th 
T-multiplet can be expressed in the U(5) basis as 

\p;N,r,n A ,L,M) = £ ^ \[N]{n d )(r)n A LM) (64) 

with rid — T i T + 2, . . . , iV — 1 or N. Similarly, the same wave functions can be expressed in the 0(6) 
basis as 

\v;N,T,n A ,L,M) = ^ r,^ \ [N] (a) (r)n A LM) (65) 

with a = N, N — 2, . . . , (r + 1) or r. The amplitudes in Eqs. (|64|) - (!65|) depend, in general, on the actual 
values of e, a, A as well as on N. 

The 0(5) symmetry of (1631) can be further used to derive special properties of electromagnetic 
transitions. The Il( 2 ) part of the E2 operator, Eq. ([2]), is a r = 1 tensor under 0(5) and connects states 
with At = ±1. The part is a r = 2 tensor under 0(5) and connects states with At = 0, ±2. The 
B(E2) values of transitions from a state with r', n' A , V in the v' multiplet to a state with r, n A , L in the 
v multiplet can be written in the form 



B(E2; u', t' , n' A , L' — > u, r, nA, L) 

t' 

77, a L' 2 



e B T^{y,v ,Tj 



n A L 



t' = t ± 1 , (66a) 



2 

' X t' = r , r ± 2 . (66b) 



n A L 



Each B(E2) is a product of two factors. The first factor which is determined by the Hamiltonian 
depends on iV and on u, t of the initial and final states. The second factor is the 0(5) D 0(3) isoscalar 
factor (ISF). It is the same for every Hamiltonian ( 16"3"|) . is completely determined by 0(5) symmetry 
and depends only on the r, n A , L of the intial and final states. Analytic expressions for some of these 
ISF are available [HE1H6], as well as a computer code for their evaluation [4"T] . 

The factorization observed in Eq. (166]) is a manifestation of the Wigner-Eckart theorem for the 0(5) 
group and has important implications with respect to E2 rates. First, consider transitions between 
states in any u, r-multiplet to states in any z/V-multiplet (including v = v'\ A direct consequence 
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of (IST)]) is that B(E2) ratios of such transitions are equal to the ratios of ISFs, hence, are independent of 
v, v' throughout the 0(6)-U(5) region. Such ratios are also independent of N and the actual parameters 
in (|63|) . i.e., they should be the same in all 0(5) nuclei considered. Second, consider E2 transitions 
between a state with given r, n A , L and a state with given t', n' A , V either with v — v' or v ^ u', and 
even in different nuclei. In B(E2) ratios of such transitions the ISFs in Eq. ( 16 6 p cancel, and the ratio 
of the other factors [e.g., e^T^iy, v' , t) j e^T N i{y" , v'" , r")} is independent of n A , L, n' A , V . Thus, such 
ratios should be the same for all transitions between the states of the r, r'-multiplets considered. 

The very definite statements made above about E2 transitions follow directly from 0(5) symmetry. 
They hold exactly for any values of N and of parameters in the Hamiltonian ( 163|) and, in particular, at 
the 0(6) and U(5) limits. They played an instrumental role in identifying empirical signatures of 0(5) 
symmetry which are common to all IBM Hamiltonians (16 3 p in the 0(6)-U(5) region and which should 
be clearly distinguished from features, e.g., absolute B(E2) values, that can yield crucial evidence for 
0(6) or U(5) symmetries or for deviations from these limits [18H50] . 

Hamiltonians with 0(5)-PDS of type II have been used extensively for studying transitional nuclei 
in the Ru-Pd [5T] and Xe-Ba regions, whose structure varies from spherical [U(5)] to 7-unstable 
deformed [0(6)]. Such 0(5)-PDS is also relevant to the coexistence of normal and intruder levels in 
112 Cd [53]. The particular Hamiltonian ifo(5) of Eq. fl63|) . with 0(5)-PDS of type II, has also selected 
U(5) basis states as eigenstates and hence exhibits also U(5)-PDS of type I. This follows from the fact 
that APqPq, which is the only term in Ho(5) that breaks the U(5) DS, involves the operator Pq = d-d—s 2 . 
The latter operator annihilates the U(5) basis states, | [N], n d = r = N,L) and | [N], n d = r = N — 1, L) , 
for reasons given after Eq. f l2l?|) . 

3.2 0(6) PDS (type II) 

An alternative situation where PDS of type II can occur is when the entire eigenspectrum retains only 
some of the symmetries Gi in a given dynamical symmetry chain fl60|) . Such a scenario was considered 
in [H] in relation to the 0(6) chain 

U(6) D 0(6) D 0(5) D 0(3) 
I I I I . (67) 

[N] (0,(7,0) (r,0) n A L 

The following Hamiltonian, with 0(6)-PDS of type II, has been proposed 

Hi = k p£p + k 2 (n< 2 ) x n( 2 )) (2) • n< 2 ) . (68) 

The k term is the 0(6) pairing term of Eq. (I5T?|) . It is diagonal in the dynamical-symmetry basis 
\[N](a)(r)n A LM) of Eq. flBT|) with eigenvalues k (N — a)(N + a + 4). The n 2 term is constructed only 
from the 0(6) generator, = d)s + s^d, which is not a generator of 0(5). Therefore, it cannot connect 
states in different 0(6) irreps but can induce 0(5) mixing subject to At = ±1,±3. Consequently, all 
eigenstates of Hi have good 0(6) quantum number a but do not possess 0(5) symmetry r. These are 
the necessary ingredients of an 0(6) PDS of type II associated with the chain of Eq. (16 7p . 

As shown in Fig. HJ a typical spectra of Hi displays rotational bands of an axially-deformed nucleus. 
All bands of Hi are pure with respect to 0(6). This is demonstrated in the left panel of Fig. [5] for the 
K = Oi, 2i, 23 bands which have a = N, and for the K = O2 band which has a = N — 2. In this case, 
the diagonal K -term in Eq. (1681) simply shifts each band as a whole in accord with its a assignment. 
On the other hand, the K 2 -term in Eq. (|6"8"|) is an 0(5) tensor with (r, 0) = (3,0) and, therefore, all 
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eigenstates of H± are mixed with respect to 0(5). This mixing is demonstrated in the right panel of 
Fig. [5] for the L = 0, 2 members of the ground band. 

A key element in the above procedure for constructing Hamiltonians with 0(6)-PDS of type II, is 
the tensorial character of the generators contained in 0(6) but not in 0(5) [44] . In the present case, the 
tensor character of the operator IT 2 ) under 0(5) is (r, 0) = (1,0) and under 0(3), L = 2. A quadratic 
interaction II (2) • Il (2) corresponds to the 0(5) multiplication (1, 0) x (1, 0) = (2, 0) © (1, 1) © (0, 0). Since 
only the irrep (0,0) contains an L = 0, it follows that the quadratic terms must be an 0(5) scalar. 
Indeed, from Table [16] of the Appendix, we find II^ • = Co{6) ~ C*o(5) • in the next cubic order, 
the interaction (Il( 2 ) x lT 2 ))( 2 ) • Il( 2 ) corresponds to (1, 0) x (1, 0) x (1, 0); 0(5) multiplication show that 
there is only one 0(3) scalar and it has 0(5) character (3,0). Consequently, (Il( 2 ) x lT 2 ))( 2 ) • II^ 2 ^ is an 
example of a cr-conserving, r-violating interaction; it mixes (r, 0) with (r ± 1,0) and (r ± 3,0). This 
discussion highlights the fact that the existence of Hamiltonians with PDS of type II, constructed in 
this manner, may necessitate higher-order terms. 

A similar procedure can been applied to the U(6) D U(5) D 0(5) D 0(3) chain of the IBM [44J. 
The generators contained in U(5) but not in 0(5) are = (d)d)^ with L = 0,2,4. is a scalar 
in 0(5) and hence does not mix 0(5) irreps. The operators and Ujy , on the other hand, have 0(5) 
character (2,0). 0(3)-scalar interactions obtained from quadratic combinations of such tensors involve 
terms of the U(5) DS Hamiltonian, Eq. (fTBl . hence do not induce 0(5) mixing among symmetric, (r, 0) 
irreps. On the other hand, cubic 0(3)-scalar combinations of llj? and uff* can lead to two independent 
d-boson interaction terms that can induce 0(5) mixing but conserve the U(5) quantum number, rid- By 
definition, such n^- conserving but r-violating cubic terms exemplify a U(5) PDS of type II. 



4 PDS type III 

PDS of type III combines properties of both PDS of type I and II. Such a generalized PDS [51] has a 
hybrid character, for which part of the states of the system under study preserve part of the dynam- 
ical symmetry. In relation to the dynamical symmetry chain of Eq. ([6]), Gdyn D G D ■ ■ ■ D G sym , 
with associated basis, |[/ijv](S)A), this can be accomplished by relaxing the condition of Eq. ([8]), 
T[h n }(a)\\ [^iv](£o)A) — 0, so that it holds only for selected states A contained in a given irrep (E ) 
of G and/or selected (combinations of) components A of the tensor Tj^j^A- Under such circumstances, 
let G' ^ G sym be a subalgebra of G in the aforementioned chain, G D G' . In general, the Hamilto- 
nians, constructed from these tensors, in the manner shown in Eq. (Q, are not invariant under G nor 
G'. Nevertheless, they do posses the subset of solvable states, |[/jjv](£ )A), with good G-symmetry, 
(S ), while other states are mixed. At the same time, the symmetry associated with the subalgebra 
G', is broken in all states (including the solvable ones). Thus, part of the eigenstates preserve part of 
the symmetry. These are precisely the requirements of PDS of type III. In what follows we explicitly 
construct Hamiltonians with such properties within the IBM framework. 



4.1 0(6) PDS (type III) 

PDS of type III associated with the 0(6) chain, Eq. (|67t) . can be realized in terms of Hamiltonians which 
have a subset of solvable states with good 0(6) symmetry but broken 0(5) symmetry. Hamiltonians 
with such property can be constructed [51] by means of the following boson-pair operators with angular 



24 



12 


— 




9— 6 




— 


8 4- 


— 


— 


io- 7 " 2 = 

6— a=N-2 




mixed a 


~~ 0= 
K=0 2 


5 — 
4 


— 




8-3- 


~~ 2 

mixed a 


— 2— 
K=2, 


6 






4 


2 


EXP 


2 


a=N 




a=N (solvable) 


162- 

o °y 

K=0 1 



Figure 4: Experimental spectra (EXP) of 162 Dy compared with calculated spectra of Hi + C\ ($0(3), Eq. ([68]) . 
with 0(6)-PDS type II and of H 2 + C 2 C ( 3) , Eq. ([75]), with 0(6)-PDS of type III. The parameters (in keV) 
are kq = 8, k 2 = 1.364, C\ = 8 and ho = 28.5, h 2 = 6.3, C 2 = 13.45, and boson number N = 15. Adapted 
from [54]. 



momentum L = 0, 2 

Pi = V2S/ + Vl{SS)f . (69b) 

From Table[5]one sees that the P pair is an 0(6) tensor with (cr = 0, r = 0, L = 0), while P 2 involves 
a combination of tensors with (a — 2, r = 1, L — 2) and (cr = 2, r = 2, L = 2) 

Pj = 2\/3-Bj 2 ]( ^ ) 0;00 , (70a) 

= v/ 2-^[2](2>(l)0;2/i + v/ ^-^[2](2)(2)0;2 M • (^0b) 



These operators satisfy 



P |c; N) = (71a) 
*V|c;iV) = (71b) 



where 

\c;N) = (N\r 1/2 (bl) N \0) , 6t = ( 4 + s t)/V2. (72) 

The state |c; AT) is obtained by substituting the 0(6) deformation, (3 = 1, as well as 7 = in the 
coherent state of Eq. i^, |c; N) = \ j5 = 1, 7 = 0; iV). It has good 0(6) character, (a) = (N), and serves 
as an intrinsic state for a prolate-deformed ground band. Rotational members of the band with a = N 
and even values of L, are obtained by angular momentum projection, \[N](N)LM) oc Vlm\c; N). The 
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Figure 5: Left: 0(6) decomposition of wave functions of states in the bands K = 0i, 2i, O2, (L = K + ), and 
= 23, (L = 3 + ), for H\ ((68D with 0(6)-PDS type II (upper portion) and H 2 ^ with 0(6)-PDS type III 
(lower portion). Right: 0(5) decomposition of wave functions of L = 0, 2 states in the a = N ground bands 
{K = 0i) of H\ (upper portion) and H2 (lower portion). Adapted from |54j . 

projection operator, Vlm, involves an 0(3) rotation which commutes with P Q and transforms among 
its various components. Consequently, P and P 2m annihilate also the projected states 

P \[N}(N}LM) = 0, (73a) 
P 2lI \[N](N}LM) = 0, L = 0,2,4,...,2iV. (73b) 

It should be noted that P\^ and P 2m , Eq. (I70bl) . span only part of the a = 2 irrep. Consequently, the 
projected states \ [N](N)LM) ofEq. ( 1731 span only part of the 0(6) irrep, (a) = (N). The corresponding 
wave functions contain a mixture of components with different 0(5) symmetry r, and their expansion 
in the 0(6) basis \[N](a)(r)n A LM) reads 

I [N] (N)LM) = < L) £ a%£> I [N] (N) (r)n A LM) , (74a) 

T,n A 

= [(N-t)\(N + t + 3)\}- 1/2 /# a . (74b) 

Here A/^^ is a normalization coefficient and explicit expressions of the factors fr L n A for L = 0, 2, 4 are 
given in Table [HJ 

Following the general algorithm, a two-body Hamiltonian with 0(6) partial symmetry can now be 
constructed [54J from the boson-pairs operators of Eq. ( |69l) as 

H 2 = h Q P^P + h 2 P}-P 2 ■ (75) 

The ho term is the 0(6)-scalar interaction of Eq. (|59|) . The multipole form of the h 2 term involves the 
Casimir operators of 0(5) and 0(3) which are diagonal in a and r, terms involving fid which is a scalar 
under 0(5) but can connect states differing by Act = 0, ±2 and a n^ 2 ^ • term which induces both 
0(6) and 0(5) mixing subject to Act = 0, ±2 and At = ±1,±3. Although H 2 is not an 0(6)-scalar, 
relations (!7ll and (1731 ensure that it has an exactly solvable ground band with good 0(6) symmetry, 
(a) = (N), but broken 0(5) symmetry. The Casimir operator of 0(3) can be added to H 2 to obtain 

Hpds = H 2 + C C (3) ■ (76) 
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Table 6: Calculated and observed (Exp) B(E2) values (in 10 2 e 2 b 2 ) for g — >• g and 7 — s- g transitions in 162 Dy. 
The parameters of the E2 operator, Eq. ©, are e B = 0.138 [0.127] eb and x = -0.235 [-0.557] for Hi (JBH} 
[J3" 2 CSW- The Hamiltonian H 1 (H 2 ) has 0(6)-PDS of type II (type III). Adapted from |5l]. 
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The solvable states of Hp^s form an axially-deformed ground band 

\[N](N)LM) L = 0,2,4, ... , 2N (77a) 
E PDS = CL(L + 1) . (77b) 

Thus, Hpds dZH]) has a subset of solvable states with good 0(6) symmetry, which is not preserved by 
other states. All eigenstates of Hp^s break the 0(5) symmetry but preserve the 0(3) symmetry. These 
are precisely the required features of 0(6)-PDS of type III. 

The calculated spectra of H 2 (1751) and H x (168]) . supplemented with an 0(3) term, are compared 
with the experimental spectrum of 162 Dy in Fig. |H The spectra display rotational bands of an axially- 
deformed nucleus, in particular, a ground band (K = 0i) and excited K — 2i and K = 2 bands. The 
0(6) and 0(5) decomposition of selected bands are shown in Fig. [5j For Hi, characteristic features of 
the results were discussed in Subsection 13.21 For H 2 , the solvable K = 0i ground band has a = N and 
all eigenstates are mixed with respect to 0(5). However, in contrast to H x , excited bands of H 2 can 
have components with different 0(6) character. For example, the K = 2 band of H 2 has components 
with a = N (85.50%), a = N - 2 (14.45%), and a = N - 4 (0.05%). These cr-admixtures can, in 
turn, be interpreted in terms of multi-phonon excitations. Specifically, the K = 2 band is composed of 
36.29% (5, 63.68% 7x =0 , and 0.03% (3 2 modes, i.e., it is dominantly a double-gamma phonon excitation 
with significant single-/? phonon admixture. The K = 2\ band has only a small 0(6) impurity and 
is an almost pure single-gamma phonon band. The results of Fig. [5] illustrate that H 2 (1751) possesses 
0(6)-PDS of type III which is distinct from the 0(6)-PDS of type II exhibited by Hi (|5gjl . 

In Table [6] the experimental B(E2) values for E2 transitions in 162 Dy, are compared with PDS 
calculations. The B(E2) values predicted by Hi (1681) and H 2 (1751) for K = 0i — > K — 0i and K = 
2i — > K — Oi transitions are very similar and agree well with the measured values. On the other 
hand, their predictions for interband transitions from the K = 2 band are very different [51]. For 
Hi, the K — O2 — > K — 0i and K = 2 — » K = 2i transitions are comparable and weaker than 
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Table 7: Calculated and observed (Exp) [55] B(E2) values (in e 2 b 2 ) for transitions from the K = O2 band in 
162 Dy. The calculations involve the Hamiltonian H2 (|75p with 0(6)-PDS of type III and the CQF Hamiltonian 
with broken 0(6) symmetry. Adapted from [55] • 
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K — 2i — > K — 0i. In contrast, for i^ 2? ^ = 2 — > K — 2j and If = 2i — >■ K — 0i transitions are 
comparable and stronger than K — 2 ~ * K — Ox. The results of a recent detailed measurement [53] 
of 162 Dy, shown in Table [7J indicate that characteristic features of the K = 2 band in this nucleus 
are reproduced by both H2 with 0(6)-PDS of type III, and the CQF Hamiltonian with broken 0(6) 
symmetry, but refinements are necessary. 



5 Partial Solvability 



The PDS of type I and III, discussed so far, involve subsets of solvable states with good symmetry 
character, with respect to algebras in a given dynamical symmetry chain. A further extension of 
this concept is possible, for which the selected solvable states are not associated with any underlying 
symmetry. Such a situation can be referred to as partial solvability. In the PDS examples considered 
within the IBM framework, the solvable states were obtained by choosing specific deformations and 
projecting from an intrinsic state, Eq. (j5]), representing the ground band 



ft,Ti N ) oc /3cos74 + /3sin7(4 + dL 2 )/v / 2 + 



N 



10) 



(7f 



Specifically, for SU(3)-PDS of type I, the solvable ground band was associated with deformations (ft = 
\/2,7 = 0), while for 0(6)-PDS of type III, it was associated with (ft = 1,7 = 0). More generally, a 
natural candidate for a solvable ground band would be the set states of good 0(3) symmetry L, projected 
from the prolate-deformed intrinsic state, | ft, 7 = 0; N), with arbitrary deformation ft > |56j 



ft; N, LM) oc 



r!v L) (/3) 



rSv L) (/3) 



'1/2 



V LM \ftn = 0;N) 

^\ f Q dx [l+ft 2 P2(x)] N P L (x) . 



0,2,4, 



2N 



(79) 



Here Pl(x) is a Legendre polynomial with L even and rffl(ft) is a normalization factor. In general, these 
L-projected states do not have good symmetry properties with respect to any of the IBM dynamical 
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Table 8: The factors /j# A , Eq. ([HID, f or the states \P;N,LM), Eq. jSD]), with L = 0,2,4. The label n A is not 
required, since these L-states are multiplicity-free. 
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symmetry chains ([2]). Their wave functions have the following expansion in the U(5) basis 

\/3;N,LM) = i[ 1 + (- 1 ) rtd " T ] Zl N JL\l N }(n d )(T)n A LM) , (80) 



2 



where (r, wa) take the values compatible with the 0(5) D 0(3) reduction and the n d summation covers 
the range t < n d < N . The coefficients £n^i A are of the form [5j 



J [(iV - n d )\(n d - T)\\(n d + r + 3)!!J ' 

Explicit expressions |57j for some of the factors fryn A are given in Table [HJ 

The construction of partially-solvable Hamiltonians which have the set of states ( 179]) as eigenstates, 
can be accomplished (58] by means of the following boson-pair operators with angular momentum 
L = 0, 2 

P f (/3 ) = S ■ S - /3 V) 2 , (82a) 
Plfa) = /9oV2dt S t + V7 . (82b) 



-1/2 f3 n d 



These operators satisfy 



or equivalently, 



Wo) I A), 7 = 0; iV> = 0, (83a) 
i\,(A))|A),7 = 0;A0 = 0, (83b) 



P (f3 )\f3 -N,LM) = 0, 

P 2 M\/3 ;N,LM) = 0. (84) 



The following Hamiltonian [58Tf60] 

H(h ,h 2 ,(3 ) = h PiWo)PoWo)+h 2 P}{(3a)- P 2 (A>) , (85) 

has a solvable zero-energy prolate-deformed ground band, composed of the states in Eq. ( 1791) . The 
Casimir operator of 0(3) can be added to it to form a partially solvable (PSolv) Hamiltonian 

Hpsoiv = H(h ,h 2 ,(3o) + CC 0{ 3) ■ (86) 
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The solvable states and energies are 

\(3 ;N,LM) L = 0,2,4,..., 2N (87a) 
E PSolv = CL(L + 1) . (87b) 

Since the wave functions of these states are known, it is possible to obtain closed form expressions 
for related observables. For example, for the E2 operator of Eq. ([3]), the B(E2) values for transitions 
between members of the solvable ground band read [56ll6Tj 

> 1 ri v L l 1 (/3) + a 2 r^ 1 2) (/3)] 2 



B(E2;L + 2^ L) = e 2 B (L + 2, 0; 2, 0|L, 0) 2 /3 2 



r^(/3)r^ +2) (^) 



a 1 = l-pJl X L/(2L + 3) , a 2 = l-/3y^ X (£ + 3)/(2L + 3) , (88) 

where the symbol (...|..) denotes an 0(3) Clebsch Gordan coefficient. 

The Hamiltonian Hpsoiv of Eq. ( 186]) is partially solvable for any value of (3q > 0. For [3q = a/2, it 
reduces to the Hamiltonian of Eq. (T4"0"]) with SU(3)-PDS of type I. In this case, the solvable states span 
the SU(3) irrep (2N, 0) and the normalization factor in Eq. (179"]) is given by 

r^ ) (/3 = v / 2) = (2iV -L)!!(2iV + ) L + l)!!iV! ' (89) 

Relation (|80|) then provides transformation brackets between these SU(3) states and the U(5) basis 
and Eq. (188]) reduces to a well-known expression for E2 transitions among states in the SU(3) ground 
band [151 EE]- When (3 = 1, the Hamiltonian Hpsoiv (|86l) coincides with the Hamiltonian of Eq. ( 1761) 
with 0(6)-PDS of type III. 

When h 2 = 0, the Hamiltonian of Eq. (|85|) takes the form 

#(/*>, A>) = hPiWPoW ■ (9°) 

Both P^(f3o) and Po(A)), Eq. (182aj) . are 0(5)-scalars. Futhermore, Po(/?o) annihilates the intrinsic state, 
Eq. (178"]) . with (3 = (3o and arbitrary 7 

p o (A))|A),7;A0 = 0. (91) 

Equivalent ly, 

P (P )\P -N,r : n A ,LM) = 0, (92) 

where the indicated states, with good r and L quantum numbers, are obtained by 0(5) projection from 
the deformed intrinsic state | {3 , 7; N) (175]) 



1/3, 7; AO 

&w = E 5 [1 + (JV -.„,„ (ra / r)!! 3 )„ • < 93 > 

Here F^\/3) is a normalization factor and the rid summation covers the range t < rid < N. The 
corresponding wave functions have the following expansion in the U(5) basis [62] 

\P;N,r,n A ,LM) = £ \ [l + (-1)^] [N]{n d )(r)n A LM) , (94a) 

) r (t) 1 _1 ^ 2 /3 nd 

^ = [ F »W)\ [(iV-n d )!(n d -r)!!(n d + r + 3)!!]V2 " (94b) 

30 



The Hamiltonian (l9"0"j) mixes the U(5) and 0(6) chains but preserves the common 0(5) subalgebra. 
This is explicitly seen from its multipole form 



H(h ,(3 ) = Mo [-C o{fi )+5N + (3 2 N(N-l) 

+h (l - ft) Un d + 2f3 2 (N - l)n d + (1 - Pl)n d (n d - 1) - Cb(5) 



(95) 



It has a solvable zero-energy 7-unstable deformed ground band, composed of the states in Eq. f|93|) . The 
Casimir operators of 0(5) and 0(3) can be added to it to form a partially solvable (PSolv) Hamiltonian 



Hpsoiv = H(h , (3 ) + B Cq(5) + C Cq(3) • 



(96) 



H psoiv (l96|) has also an 0(5)-PDS of type II in the sense discussed in Subsection 13.11 The solvable 
states and energies are 



I ft; N, t, n A , LM) , 

EpsoIv = B r(r + 3) + C L(L + 1) 



(97a) 
(97b) 



where the (r, n A , L) assignments are the same as for states in the 0(6) irrep with a = N. Closed form 
expressions can be derived for observables in these states. For example, for the general E2 operator of 
Eq. ((3D, we find [62] 



B(E2;r + l,n' A ,L' r, n A ,L) 



t + 1 o2 [FS,%(P) + F^>(/3) 



2r 



(t+1), 



T+1 1 

n' A L' 2 



T 

n A L 



(9* 



This expression is similar in form to that encountered in Eq. (I66ap . but now the factor in front of the 
0(5) isoscalar factor is explicitly known. The Hamiltonian of Eq. ( 196]) is partially solvable for any value 
of (3q > 0. For /3q = 1, it reduces to the Hamiltonian of Eq. ( 15 2 p with 0(6) dynamical symmetry. The 
solvable states ( 19 7p then span the 0(6) irrep (a) = (N) and the normalization factor (|9"3"|) becomes 



2 N+1 (N + 1)\ 
(N - t)\(N + t + 3)! ' 



(99) 



In this case, relation (IMjl corresponds to known transformation brackets between these 0(6) states and 
the U(5) basis, and one recovers from Eq. (|98|) a familiar expression for the indicated B(E2) in the 
0(6) limit of the IBM [16]. 

In addition to the states shown in Eq. ( 1921) . the operator Po(A)) (I82aj) annihilates also the following 
U(5) basis states 



P Q \[N],n d = r = N,n A ,LM) 
P \[N),n d = T = N-l,n A ,LM) 



, (100a) 
, (100b) 

for reasons explained after Eq. (T2"31) . Consequently, H PSo iv of Eq. (|9"B"1) has also a U(5)-PDS of type I, 
in the sense discussed in Subsection 12.11 The additional solvable eigenstates and energies are 

I [N] ,n d = t = N, n A , LM ) E Psolv = B N{N + 3) + C L(L + 1) , (101a) 

\[N],n d = r = N-l,n A ,LM) E Psolv = B (N - 1)(N + 2) + C L{L + 1) , (101b) 
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where the (r,n^,L) assignments are those of the 0(5) D 0(3) reduction. 

The Hamiltonian H(ho, h 2 , A>) of Eq. ( 185]) is a prototype of an intrinsic Hamiltonian which generate 
band-structure [58-60J. Its energy surface, defined as in Eq. 01]), 

(102) 



£7*00,7) = iV(iY-l)(l + /3 2 )- 2 [/ io (/3 2 -/3 2 ) 2 + 2M 2 (/3 2 -2/3o/3cos3 7 + /3 2 )] 



has a global minimum at (/3 > 0,7 = 0), corresponding to a prolate-deformed shape. H(h , h 2 , (3 ) 
is 0(3)-invariant, but has the deformed equilibrium intrinsic state, |/?o,7 = 0; N) (!7H|) . as a zero- 
energy eigenstate. The 0(3) symmetry is thus spontaneously broken. The two Goldstone modes are 
associated with rotations about directions perpendicular to the symmetry axis. The intrinsic modes 
involve the one- dimensional /3 mode and two-dimensional 7 modes of vibrations. For large N, the 
spectrum of H(ho, h 2 , 0o) (155]) is harmonic, involving and 7 vibrations about the deformed minimum 
with frequencies given by [581160] 

(103) 



e p = 2N(3l(2h Q + h 2 ) , e 7 = 18iV/3 2 (l + /3 2 )- 1 /i 2 . 



The importance of H(ho, h 2 , 0o) lies in the fact that the most general one- and two-body IBM Hamil- 
tonian with equilibrium deformations (/?o > 0,70 = 0), can be resolved into intrinsic and collective 
parts [591160] 



Hibm = H(h ,h 2 , (3 ) + H c 



(104) 



The intrinsic part is the partially-solvable Hamiltonian of Eq. ( 185]) . The collective part, H c , involves 
kinetic rotational terms which do not affect the shape of the energy surface 
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(105) 



The various Casimir operators in Eq. ( I105j) are defined in the Appendix. The L-projected states, 
I (3;N,LM), of Eq. (179]) can now be used to construct an £-projected energy surface, E^\/3) = 
{Pi N, LM\H IBM \(3] N, LM), for the IBM Hamiltonian (IT04]) 



E[ N \(3) = h (/3 2 - (3lf + 2/12 09 - A)) 2 Eg? + c 3 \L(L + 1 
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6< y 
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(106) 



Here d[ N ^, S 2 N l , an d £22? denote the expectation values in the states \(3; N, LM) of n^, n s (n s — 1), 
ndifid — 1) and h s hd respectively. All these quantities are expressed in terms of the expectation value 
of n s , denoted by sg?. Specifically, = N- S™, Sg? = fig?^ E$ = (tf - l)fig? - sg?, 
flji = iV(iV — 1) — 2(iV — 1)5*1^ + S 2L . The quantity Sf^ itself is determined by the normalization 
factors of Eq. ([79]) 

S[ N L ] = (P;N,LM\n s \P;N,LM) = r^tf) /T%\/3) . (107) 
It also satisfies the following recursion relation [56] 



°1,L 
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For h 2 = 0, the intrinsic part of the IBM Hamiltonian in Eq. (I104p reduces to the Hamiltonian 
H(ho,flo) of Eq. (1901 . which is 0(5)-invariant. The unprojected energy surface (I102p is now inde- 
pendent of 7 and the equilibrium shape is deformed ((3q > 0) and 7-unstable. The 0(5) symmetry 
is spontaneously broken in the intrinsic state, \^ ,j;N) (!7H|) . which is a zero-energy eigenstate of 
H(h , f3 ). As a result, the 7 and three rotational modes are Goldstone modes, and only the (3 vibration 
in Eq. (11031) survives as a genuine mode [SHIED]- In this case, H IBM {h 2 = 0) in Eq. (I104p preserves the 
0(5) symmetry, r, and the 0(5)-projected states of Eq. fl93l) can be used to construct its r-projected 
energy surface, E ( T N L \f3) = (/?; N, r, n A , LM\H IBM (h 2 = 0) |/5; N, r, n A , LM), 



E%!W) = h (F-ftyS£> + C3 L(L + 1) — QD\ 



2\2 o(N) 



+ C 5 



r(r + 3) - AD^J 
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+ E 



(109) 



Here d[ N J and S 2 l' denote the expectation values in the states \(3; N,r,n A , LM) of and h s (h s — 1) 
respectively. All these quantities are expressed in terms of the expectation value of h s , denoted by 
. Specifically, D[\} = N — S^, = S^) S^T 1 ' . The quantity itself is determined by the 
normalization factors of Eq. (I9"3"|) 
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S$ = (/3;iV,r,n A ,LM|n s |/3;iV,r,n A ,LM) = F^P) / (?) . 
It also satisfies the following recursion relation 

y (N) _ (N-t)(N + t + 3) 
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6 PDS and Quantum Phase Transitions 

Symmetry plays a profound role in quantum phase transitions (QPT). The latter occur at zero tem- 
perature as a function of a coupling constant in the Hamiltonian. Such ground-state energy phase 
transitions [63] are a pervasive phenomenon observed in many branches of physics, and are realized 
empirically in nuclei as transitions between different shapes. QPTs occur as a result of a competi- 
tion between terms in the Hamiltonian with different symmetry character, which lead to considerable 
mixing in the eigenfunctions, especially at the critical-point where the structure changes most rapidly. 
An interesting question to address is whether there are any symmetries (or traces of) still present at 
the critical points of QPT. As shown below, unexpectedly, partial dynamical symmetries (PDS) can 
survive at the critical point in spite of the strong mixing [61J. The feasibility of such persisting symme- 
tries gains support from the recently proposed [61] and empirically confirmed [62] analytic descriptions 
of critical-point nuclei, and the emergence of quasi-dynamical symmetries [27] in the vicinity of such 
critical-points. 

A convenient framework to study symmetry-aspects of QPT in nuclei is the IBM [8] , whose dynamical 
symmetries (T5]) correspond to possible phases of the system. The starting point is the energy surface of 
the Hamiltonian, Eq. which for one- and two- body interactions has the form 

E N (P,-y) = £ + iV(iV-l)/(/3,7) , 
/(ft 7) = (l + /3 2 )~ 2 /3 2 [a-6/3cos3 7 + c/3 2 ] . (112) 

The coefficients E ,a,b,c involve particular linear combinations of the Hamiltonian's parameters [60J. 
Phase transitions can be studied by IBM Hamiltonians of the form, H(a) = (1 — a) Hi + a H 2 , involving 
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Figure 6: Energy surfaces at the critical points, Eq. (I113j) . (a) First-order transition. The position and height 
of the barrier are ft = ft + = (-1 + y / T+ft^ )/fto and h = /(/?+, 7 = 0) = (-1 + y/l + ft' 2 . ) 2 /4 respectively, 
(b) Second-order transition. In this case f(ft, 7) is independent of 7. Asymptotically, f(ft — > 00,7) = 1. 



terms from different dynamical symmetry chains [20]. The nature of the phase transition is governed 
by the topology of the corresponding surface fll!2j) . which serves as a Landau's potential with the 
equilibrium deformations as order parameters. The conditions on the parameters and resulting surfaces 
at the critical-points of first- and second-order transitions are given by 

1 st order b 2 = 4ac, a > 0, b > f(ft, 7 = 0) = c(l + ft 2 )- 2 ft 2 (ft - f3 ) 2 , (113a) 
2 nd order a = 0, 6 = 0, c>0 /(/?, 7) = c(l + /3 2 )" 2 /3 4 . (113b) 

As shown in Fig. [6j the first-order critical-surface has degenerate spherical and deformed minima 
at ft = and (ft = (3q > 0,7 = 0), where ft = 2a/b. The position (ft + ) and height (h) of the 
barrier are indicated in the caption. The second-order critical-surface is independent of 7 and is flat 
bottomed (~ /3 4 ) for small ft. The conditions on a, b, c in Eq. (11131) fix the critical value of the control 
parameter (a = a c ) which, in turn, determines the critical-point Hamiltonian, H cri = H(a = a c ). 
IBM Hamiltonians of this type have been used extensively for studying shape-phase transitions in 
nuclei [2lJ t[2Tl[o^[oTl[oT| - lo r 7] . We now show that a large class of such critical-point Hamiltonians exhibit 
PDS [ST]. 

The spherical to deformed 7-unstable shape-phase transition is modeled in the IBM by the Hamil- 
tonian 

H m = eh d + A[d^S-(s^) 2 ] [H.c] 

e = A(N-1)A. (114) 

The A-term is the 0(6) pairing term of Eq. ( 1591) . H cr i satisfies condition ( 1113bl) with c = 4A, hence 
qualifies as a second-order critical Hamiltonian. It involves a particular combination of the U(5) and 
0(6) Casimir operators, hence is recognized to be a special case of the Hamiltonian Ho(5) of Eq. (!63l) 
with 0(5)-PDS of type II. In fact, since 0(5) is a good symmetry common to both the U(5) and 0(6) 
chains (J62l) . the 0(5) PDS is valid throughout the U(5)-0(6) transition region. As mentioned at the 
end of Subsection 13.11 i/o(5) an d, therefore, H cr i fl!14j) . has also U(5)-PDS of type I, with the following 
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64320 023456789 10 L n d n d 

Figure 7: Left: spectrum of a first-order critical Hamiltonian Hcri(fio), Eq. (|116p . with h 2 = 0.05, (3q = 1.3 
and N = 10. The solvable eigenstates are the deformed states, Eq. (|117p . forming a zero-energy K = 0i 
ground band and the spherical states, L = O2, 3i, Eq. f j 1 1 9 [) . with good U(5) symmetry. Right: U(5) {n d ) 
decomposition for selected eigenstates of H cr i(/3o). Adapted from [56J. 



solvable U(5) basis states 

\[N],n d = t = N,L) E = eN, (115a) 

|[JV],7i ( , = r = JV-l,L) E = e(N-l), (115b) 

where L takes the values compatible with the 0(5) D 0(3) reduction. 

The dynamics at the critical point of a spherical to prolate-deformed shape-phase transition can be 
modeled in the IBM by the following Hamiltonian [56] 

HcrM = h 2 P^f3 )-P 2 (f3 ) , (116) 

where P 2 \{Po) is the L = 2 boson-pair of Eq. (I82bl) and h 2 , (3q > 0. The corresponding surface in 
Eq. (11121) has coefficients a = 2h 2 f3^,b = 4h 2 (3o,c = 2h 2 , which satisfy condition flll3a[) . This qualifies 
H C ri{Po) as a first-order critical Hamiltonian whose potential accommodates two degenerate minima 
at j3 = and (/?, 7) = ((3 ,0). H cri ((3 ) is recognized to be a special case of the partially-solvable 
Hamiltonian, Hpsoiv of Eq. (156]) . As such, it has a solvable prolate-deformed ground band, composed 
of the states of Eq. (187)1 

\/3 ;N,L) E = L = 0,2,4, ... , 2N . (117) 
On the other hand, the following multipole form of H cr i((3o) 

H cr i(f3o) = 

h 2 [2(f3 2 N - 2)h d + 2(1 - Pl)h 2 d + 2(70(5) - ^0(3) + y/UfaU® ■ f/ {2) ] (118) 

identifies it as the Hamiltonian of Eq. ( IT8|) with U(5)-PDS of type I. As such, it has also the solvable 
spherical eigenstates of Eq. (I2ip . with good U(5) symmetry 

\N,n d = r = L = 0) E = (119a) 
\N, n d = t = L = 3) E = 6h 2 [/3$(N - 3) + 5] . (119b) 
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Figure 8: Left: spectrum of i? cri (/3 = \/2), Eq. (fTIS) . with /t 2 = 0.05 and iV = 10. L(K = 0i) and 
L(if = 2i) are the solvable SU(3) states of Eq. (fT20al) and Eq. (fT20E]l with fc = 1, respectively. L = 2 , 3i are 
the solvable U(5) states of Eq. (|119p . Right: U(5) (rid) and SU(3) [(A,//)] decomposition for selected eigenstates 
of H cri {/3 = a/2). Adapted from [SI]. 



The spectrum of H cr t(f3o) (I116P and the U(5) (rid) decomposition of selected eigenstates is shown in 
Fig. [7] The spectrum displays a coexistence of spherical states (some of which solvable with good 
U(5) symmetry) and deformed states (some of which solvable), signaling a first-order transition. The 
remaining non-solvable states in the spectrum are either predominantly spherical (with characteristic 
dominance of single components) or deformed states (with a broad rid distribution) arranged in 
several excited bands [56] . 

The critical Hamiltonian of Eq. (I116P with /3 = v^2 is a special case of the Hamiltonian of Eq. fpfO]) . 
shown to have SU(3)-PDS of type I. As such, it has a subset of solvable states, Eqs. ()4Tl)-fj42l). which 
are members of the ground g(K = 0) and / ~f k (K = 2k) bands, with good SU(3) symmetry, (A,/i) = 
(2N - 4k, 2k) 

\N, (2N,0)K = 0,L) E = L = 0,2,4, ... , 2N (120a) 
\N, (2N - 4k, 2k)K = 2k, L) E = h 2 Qk (2N -2k + 1) 

L = K,K+l,...,{2N -2k) k > . (120b) 

In addition, H cri ((3 = a/2) has the spherical states of Eq. ( 11 1911 . with good U(5) symmetry, as eigen- 
states. The spherical L = state, Eq. f)119ap . is exactly degenerate with the SU(3) ground band, 
Eq. (I120al) . and the spherical L = 3 state, Eq. (I119bl) . is degenerate with the SU(3) 7-band, Eq. (112Qb|) 
with k = 1. The remaining levels of H cr i((3o = a/2), shown in Fig. [8] are calculated numerically and 
their wave functions are spread over many U(5) and SU(3) irreps. This situation, where some states are 
solvable with good U(5) symmetry, some are solvable with good SU(3) symmetry and all other states 
are mixed with respect to both U(5) and SU(3), defines a U(5) PDS of type I coexisting with a SU(3) 
PDS of type I. 

The critical Hamiltonian of Eq. (1 1 1 6 j) with /3 = 1 is a special case of the Hamiltonian of Eq. (!76|) . 
shown to have 0(6)-PDS of type III. As such, it has a subset of solvable states Eq. ( JT71) . which are 
members of a prolate-deformed ground band, with good 0(6) symmetry, (a) = (N), but broken 0(5) 
symmetry 

\N,a = N,L) E = L = 0,2,4, ... , 2N . (121) 
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Figure 9: Upper panel: spectrum of H cr i{Po = 1), Eq. (|116[) . with h 2 = 0.05 and N = 10. L(K = 0i) are the solvable 
states of Eq. (|121[) with good 0(6) but broken 0(5) symmetry. L = 02,3i are the solvable U(5) states of Eq. (|119[l . 
Bottom left panel: U(5) (rid) and 0(6) (<r) decomposition for selected spherical and deformed eigenstates of H cr i(f3o = 1)- 
Bottom right panel: 0(5) (r) decomposition for the L — 0, 2 states, Eq. (|121|) . members of the ground band (K = 0i) of 
H cr i((3o — 1). Both states have 0(6) symmetry a — N. Adapted from |61j . 



In addition, H cri ((3o = 1) has the spherical states of Eq. f 1 1 1 9 [) . with good U(5) symmetry, as eigenstates. 
The remaining eigenstates of H cr i((3o = 1) shown in Fig. [9] are mixed with respect to both U(5) and 
0(6). Apart from the solvable U(5) states of Eq. f 1 11 9 [) . all eigenstates of H cr i((3o = 1) are mixed with 
respect to 0(5) [including the solvable 0(6) states of Eq. (112 ip . as shown in the bottom right panel of 
Fig. [9]. It follows that the Hamiltonian has a subset of states with good U(5) symmetry and a subset of 
states with good 0(6) but broken 0(5) symmetry, and all other states are mixed with respect to both 
U(5) and 0(6). These are precisely the required features of U(5) PDS of type I coexisting with 0(6) 
PDS of type III. 

In conclusion, the above results demonstrate the relevance of the PDS notion to critical-points of 
QPT, with phases characterized by Lie-algebraic symmetries. In the example considered, second-order 
critical Hamiltonians mix incompatible symmetries but preserve a common lower symmetry, resulting 
in a single PDS with selected quantum numbers conserved. First-order critical Hamiltonians exhibit 
distinct subsets of solvable states with good symmetries, giving rise to a coexistence of different PDS. 
The ingredients of an algebraic description of QPT is a spectrum generating algebra and an associated 
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geometric space, formulated in terms of coherent (intrinsic) states. The same ingredients are used in the 
construction of Hamiltonians with PDS. These, in accord with the present discussion, can be used as 
tools to explore the role of partial symmetries in governing the critical behaviour of dynamical systems 
undergoing QPT. 

7 PDS and Mixed Regular and Chaotic Dynamics 

Partial dynamical symmetries can play a role not only for discrete spectroscopy but also for analyzing 
statistical aspects of nonintegrable systems p8j|69] . Hamiltonians with dynamical symmetry are always 
completely integrable [70]. The Casimir invariants of the algebras in the chain provide a set of constants 
of the motion in involution. The classical motion is purely regular. A dynamical symmetry-breaking 
is connected to nonintegrability and may give rise to chaotic motion [70H72"] . Hamiltonians with PDS 
are not completely integrable, hence can exhibit stochastic behavior, nor are they completely chaotic, 
since some eigenstates preserve the symmetry exactly. Consequently, such Hamiltonians are optimally 
suitable to the study of mixed systems with coexisting regularity and chaos. 

The dynamics of a generic classical Hamiltonian system is mixed [73] ; KAM islands of regular motion 
and chaotic regions coexist in phase space. In the associated quantum system, if no separation between 
regular and irregular states is done, the statistical properties of the spectrum are usually intermediate 
between the Poisson and the Gaussian orthogonal ensemble (GOE) statistics. In a PDS of type I, the 
symmetry of the subset of solvable states is exact, yet does not arise from invariance properties of the 
Hamiltonian. This offers an important opportunity to study how the existence of partial (but exact) 
symmetries affects the dynamics of the system. If the fraction of solvable states remains finite in the 
classical limit, one might expect that a corresponding fraction of the phase space would consist of KAM 
tori and exhibit regular motion. It turns out that PDS has an even greater effect on the dynamics. It 
is strongly correlated with suppression (i.e., reduction) of chaos even though the fraction of solvable 
states approaches zero in the classical limit [HBJEH]- 

We consider the IBM Hamiltonian of Eq. (I85j) 

ff(A>) = ho P t (/3 )P (/3o) + h P 2 f (/3 ) • F 2 (A>) • (122) 

As discussed in Section [51 when 0q = y/2, the Hamiltonian f 11 2 2 j) has an SU(3)-PDS of type I. In this 
case, the solvable states are those of Eqs. (|4T]) - (l42|) . At a given spin per boson I = L/N, and to leading 
order in 1/N, the fraction / of solvable states decreases like 1/N 2 with boson number. However, at 
a given boson number N, this fraction increases with /, a feature which is valid also for finite iV [68J. 
The classical limit of HI 22 j) is obtained [73H76] through the use of coherent states parametrized by the 
six complex numbers {a s , a M ; /i = —2, . . . , 2} and taking iV — > oo. The classical Hamiltonian is then 
obtained from (II 2 2 p by the substitution s\ — > a*, a* and s, — > a s , and rescaling the parameters 
hi Nhi (i = 0, 2). Here 1/N plays the role of H. 

To study the effect of the SU(3) PDS on the dynamics, we fix the ratio h 2 /ho at a value far from the 
exact SU(3) symmetry (for which ho/h 2 = 1). We then change (3q in the range 1 < /3o < 2. Classically, 
we determine the fraction a of chaotic volume and the average largest Lyapunov exponent A. To 
analyze the quantum Hamiltonian, we study spectral and transition intensity distributions. The nearest 
neighbors level spacing distribution is fitted by a Brody distribution, P aJ (S') = AS W exp(— aS 1+UJ ), where 
A and a are determined by the conditions that Pui(S) is normalized to 1 and (S) = 1. For the Poisson 
statistics u = and for GOE u = 1, corresponding to integrable and fully chaotic classical motion [TTJ 
[78], respectively. The intensity distribution of the SU(3) E2 operator, of Eq. (H4l) . is fitted by a x 2 
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Figure 10: Classical (a, A) and quantal (oj, v) measures of chaos versus j3o for the Hamiltonian (|122p with 
^2/^-0 = 7.5. Shown are three cases with classical spins I = 0.08, 0.4, and 1. The quantal calculations (u, v) 
are done for N = 25 bosons and spins L = 2, 10, and 25, respectively. Notice that with increasing spin the 
minimum gets deeper and closer to /3o = y/2. The suppression of chaos near /?o = v2 is seen both for finite N 
through the measures to, v and in the classical limit N — > 00 through the measures a, A. Adapted from 



distribution in v degrees of freedom [79], P v {y) = l(v/2(y)y /2 /T(b>/2)]y u / 2 ~ l exp(—uy/2(y)). For the 
GOE, v — 1 and v decreases as the dynamics become more regular. 

Fig. [TOl shows the two classical measures a, A and the two quantum measures u, v for the Hamilto- 
nian (I122p as a function of (3q- The parameters of the Hamiltonian are taken to be ^2/^0 — 7.5 and the 
number of bosons is N = 25. Shown are three classical spins I = 0.08, 0.4 and 1, which correspond in 
the quantum case to L = 2, 10 and 25. All measures show a pronounced minimum which gets deeper 
and closer to /3q = y2 [where the partial SU(3) symmetry occurs] as the classical spin increases. This 
behaviour is correlated with the fraction of solvable states (at a constant N) being larger at higher I. 
We remark that the classical measures show a clear enhancement of the regular motion near /3 = v2 
even though the fraction of solvable states vanishes as 1/A^ 2 in the classical limit N — > 00. 

To confirm that the observed suppression of chaos is related to the SU(3) PDS, we employ the 
concept of an entropy [80, 81J associated with a given symmetry. To determine the SU(3) entropy, we 
expand any eigenstate \aLM) in an SU(3) basis, \aLM) = 52(a^) k c (a^)kK^' l*)KLM). Denoting by 
p^ the probability to be in the SU(3) irrep (A,/i), p^ = Y^k \ c ^{\h)k\ 2 i ^ ne SU(3) entropy of the state 
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Figure 11: The average SU(3) entropy of the eigenstates of the Hamiltonian (|122p (for hijh§ = 7.5) versus 
/3o, for three values of the spin (per boson), I = 0.08, 0.4, and 1. Left: N = 15 bosons; right: N = 25 bosons. 
Adapted from [68] . 



\aLM) is defined as S^j^ = — Ylx^Px^ ^Px/l- The entropy vanishes when the state has a good SU(3) 
symmetry. The averaged entropy (Ssu(3)) over an eigenstates is then a measure of the global SU(3) 
symmetry. This quantity is plotted in Fig. [TTJ versus (3q for N = 15 and 25 and for the same spin 
values (per boson) / as in Fig. [TUl We observe a minimum which is well correlated with the minimum 
in Fig. [10J The maximum SU(3) entropy is the logarithm of the number of allowed SU(3) irreps for 
the given iV and I. The average SU(3) entropy therefore increases with N. The depth of the minimum 
increases with A" and I though the fraction of solvable states is smaller at A^ = 25 than at N = 15 by 
a factor of about 3. The existence of an SU(3) PDS seems to have an effect of increasing the SU(3) 
symmetry of all states, not just those with an exact SU(3) symmetry [68] . 

In order to better understand the strong suppression of classical chaos induced by PDS, we consider 
a simpler model and use its PDS to infer relationships between the classical and quantum dynamics of 
a Hamiltonian in a mixed KAM regime [69]. The model is based on a U(3) spectrum generating algebra 
and its building blocks are three types of bosons a', b* , c* satisfying the usual commutation relations. 
The nine number- conserving bilinear products of creation and destruction operators comprise the U(3) 
algebra. The conservation of the total boson-number A" = h a + hf, + h c {h a = a^a with eigenvalue n a etc.) 
ensures that the model describes a system with only two independent degrees of freedom. All states 
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of the model are assigned to the totally symmetric representation [N] of U(3). One of the dynamical 
symmetries of the model is associated with the following chain of algebras 



U(3) D U(2) D U(l) (123) 

Here U(2) = SU(2) x U a f,(l) with a linear Casimir n ab = n a + h b [which is also the generator of U b(l)]. 
The generators of SU(2) are J + = tfa, J_ = alb, J z = (n b — n a )/2 and its Casimir J 2 = n ab (h ab + 2)/4. 
The subalgebra U(l) in Eq. (11231) is composed of the operator J z . A choice of Hamiltonian with a U(2) 
dynamical symmetry is 

H = u a a)a + uj b b ] b = n ab -2AJ z (124) 

where u a ^ b — 1 ± A, and A is introduced to break degeneracies. Diagonalization of this Hamiltonian 
is trivial and leads to eigenenergies E na ^ nb = u a n a + u b n b and eigenstates \n a , n b , n c ) or equivalently 
\N, J, J z ) where the label J = n ab /2 identifies the SU(2) irrep. These are states with well defined n a , 
n b and n c = N — n a — n b . To create a PDS we add the term 

H 1 = b\b ] a + b ] c + a)b + c ] b)b , (125) 

which preserves the total boson number but not the individual boson numbers, so it breaks the dynamical 
symmetry. However, states of the form \n a ,n b = 0, n c ) (or equivalently \N, J = n a /2, J z = — J) ) with 
n a = 0, 1,2, . . . N are annihilated by Hi and therefore remain eigenstates of Hq + BH\. The latter 
Hamiltonian is not an SU(2) scalar yet has a subset of (N + 1) "special" solvable states with SU(2) 
symmetry, and therefore has PDS. There is one special state per SU(2) irrep J = n a /2 (the lowest 
weight state in each case) with energy u a n a independent of the parameter B. Other eigenstates are 
mixed. Although Ho and H\ do not commute, when acting on the "special" states they satisfy 

Hq, Hi~\\n a ,n b = Q,n c ) = . (126) 

To break the PDS we introduce a third interaction 

H 2 = a ] c + c ] a + 6 f c + c ] b . (127) 

The complete Hamiltonian is then 

H = Ho + BHi + C H 2 . (128) 

For B = C = we have the full dynamical symmetry; for B ^ 0, C = we have partial dynamical 
symmetry and for C / we have neither. 

The classical Hamiltonian H c i is obtained from (11281) by replacing (a', 6' , c') by complex c-numbers 
(a*,/3*,7*) and taking iV — > oo. The latter limit is obtained by rescaling B = NB, a — » a/y/N etc. 
and considering the classical Hamiltonian per boson % = Hd/N. In the present model the latter has 
the form 

U = Uo + BUx + CUi ■ (129) 

Number conservation imposes a constraint a* a + (3* (3 + 7*7 = 1, so that the phase space is compact 
and four-dimensional with a volume 2tt 2 . The total number of quantum states is (N + 1)(N + 2)/2. 
Assigning, to leading order in N, one state per {2nh) 2 volume of phase space, we identify h = l/N, so 
that the classical limit is iV — > 00. 
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Figure 12: Classical (fi) and quantum (w) measures of chaos [denoted by (•) and (x), respectively] versus C 
for the Hamiltonian (|128h with B = 0.5. Adapted from |69| . 



In all calculations reported below [69] we take A = 0.8642 and N = 60. As a first step, we fix 
B = 0.5 and vary C. As previously done, for the classical analysis we randomly sample the phase space 
and determine the fraction fi of chaotic volume (same as a in Fig. [TO]) . For the quantum analysis we 
evaluate the energy levels, calculated the nearest neighbors level spacing distribution of the unfolded 
spectrum and fitted it to a Brody distribution. The latter is specified by the fit parameter u, mentioned 
above. As shown in Fig. [12], both of these measures indicate a suppression of chaos near C = similar 
to the results of Fig. [TO] To appreciate the strong effect of the PDS (at C = 0) on the underlying 
dynamics, it should be noted that the fraction of the solvable states \n a , n& = 0, n c ) is 2/ (N + 2), which 
approaches zero in the classical limit. To measure the extent to which each eigenstate |\?) has SU(2) 
symmetry, we define variances of = (^nf — (\&|nj|\I/) 2 (i = a,b). A state which belongs to just 
one irrep of SU(2) (with well defined J, J z ) has zero variances, while a mixed state has large variances. 
These variances have the same physical content as the entropies considered before. It is instructive 
to display the average (n a ) and variance of each state, as done in Fig. [131 SU(2) PDS is present in 
Fig. [TBTa) (B 7^ 0, C = 0), Fig. [TBT b) is a blow up of Fig [TBT a). and in Fig. [TBT c) the symmetry is 
completely broken (C ^ 0). In Figs. [TBT a) and [TBT b) we see states with zero variance. These are just 
the special N + 1 states (n& = 0) discussed before, which preserve the SU(2) symmetry. In addition, 
we see families of states with small variance and small (n&) which suggests that the presence of PDS 
increases the purity of states other than the special ones. By contrast, in Fig. [TBT c) we see no particular 
structure because of the destruction of the PDS for C ^ 0. 

Considerable insight is gained by examining the classical phase space structure in terms of action- 
angle variables a = \fT a exp(— i6 a ), ft = y/Jlexp(— i9 b ) and 7 = \fT c = a/1 — J a — J a . The 9 a = —n/2 
Poincare section is shown in Fig. 14 for energy E = 1.0. When SU(2) PDS is present (B 7^ 0, C = 0), 
we see in Figs. 14(a) and 14(b) a torus with Jj, = 0, and additional perturbed tori in its neighborhood 
(small Jb). This structure is absent when the symmetry is completely broken (C 7^ 0), as shown in 
Fig. 14(c). The features in Fig. 14 persist also at other energies. To understand them, we recall that 
for B — C — 0, the Hamiltonian (11291) is integrable and all trajectories wind around invariant tori. By 
standard torus quantization (without turning points) the actions are quantized as Jj = riih = rii/N 
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Figure 13: (Left panel). The values of (n a ) and of the variance <7& (denoted by o) of each eigenstate of the Hamilto- 
nian (|128[) . (a) B = 0.5, C = (partial dynamical symmetry), (b) a blow up of (a) with superimposed results [denoted by 
(+)] of quantum perturbation theory. The families of states with low Ob have small values of (rib). ( c ) B = 0.3, C = 0.5 
(broken symmetry). Adapted from |69) . 

Figure 14: (Right panel). Poincare sections Jb versus 9b at energy E = 1.0. (a) B = 0.5, C = (partial dynamical 
symmetry), (b) a blow up of (a) with superimposed results (dashed curves) of classical perturbation theory, (c) B = 0.3, 
C = 0.5 (broken symmetry). Adapted from [69] . 



(i = a,b). In the integrable limit, quantum states are associated with toroidal manifolds in phase space. 
In case of a partial symmetry (B ^ 0, C = 0) we are led by analogy with Eq. ( I126P to seek manifolds 
M. in phase space on which 

{% Q ,Ui) =0 (130) 

v J M 

vanishes even though the Poisson bracket is not zero everywhere. In addition, we demand that 
{{Ho, Hi}, 'Ho+'Hi}\m — (in analogy to the quantum relation [[H ,Hi] , H +Hi]\n a ,n b = 0, n c ) = 0) 
so that a trajectory starting on Ai remains on Ai. The solution to these conditions is the manifold 
Jb = (3* (3 = 0, which may be interpreted as a (degenerate) torus of the H Hamiltonian. It is also a 
stable isolated periodic orbit of Hq + BHi. Quantization of the torus with J b = proceeds exactly as 
before, so we correctly predict no change in the quantum energies associated with it. The manifold Ai 
(Jb = 0) is the direct classical analogue of the special quantum states |n& = 0). It refers, however, to a 
region of phase space of measure zero, and so cannot by itself explain the observed (global) suppression 
of chaos. However, as suggested by Fig. 14, the presence of PDS induces a quasi-regular region foliated 
by tori in the vicinity of the special torus. The dynamics on a finite measure of phase space can be 
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understood by performing a perturbative calculation in the neighbourhood of A4 

For the classical perturbation calculation we set C = in Eq. f)129p and treat B as an expansion 
parameter, assuming BT-Li in Eq. H129[) is small in the neighbourhood of the special periodic orbit. 
The second order correction reproduces well the perturbed tori on the Poincare sections as shown in 
Fig. 14(b). The variances can be calculated in quantum perturbation theory. In Fig. ITBT b) we show the 
results [denoted by (+)] of the quantum perturbation theory (to order B 5 ). We see that the first few 
families of states are reproduced. It is these states which we can recover from perturbation theory and 
whose approximate symmetry is induced by the symmetry of the special states. 

The following physical picture emerges from the foregoing analysis. Near the special orbit, there 
are KAM tori, some of which are quantized. The quantum eigenstates lie on these tori, so knowing 
the classical variance of the actions of the tori tells us the variances of the states themselves, in the 
semiclassical limit. Large variances indicate the extent to which the corresponding states fail to respect 
the symmetry. This provides a measure for a separation of regular and irregular levels, as conceived 
in [22]. In the present model, the quantum states can be grouped into three classes: i) the special 
states, which observe the symmetry; ii) the "almost special states" which are accessible by perturbation 
theory; iii) the rest of the states, which are mixed. As in [73], the frontier between regular states (sets 
(i) and (ii) ) and irregular states (set (iii) ) is not sharp. 

The above discussion illustrates the effect of PDS on the quantum and classical dynamics of a 
mixed system. At the quantum level, PDS by definition implies the existence of a "special" subset 
of states, which observe the symmetry. The PDS affects the purity of other states in the system; in 
particular, neighboring states, accessible by perturbation theory, possess approximately good symmetry. 
Analogously, at the classical level, the region of phase space near the "special" torus also has toroidal 
structure. As a consequence of having PDS, a finite region of phase space is regular and a finite fraction 
of states is approximately "special" . This clarifies the observed suppression of chaos. Based on these 
arguments and the above results, it is anticipated that the suppression of chaos will persist in higher 
dimensional systems with PDS. 

8 PDS and Higher-Order Terms 

In applications of algebraic modeling to dynamical systems, there is occasionally a need, based on 
phenomenological and/or microscopic grounds, to include higher-order terms in the Hamiltonian. For 
example, in the IBM, accommodating rigid triaxial shapes and describing large anharmonicities in 
excited bands, requires at least cubic terms in the boson Hamiltonian. From a microscopic point 
of view, many-body boson interactions in the IBM, are generated by the mapping of fermion pairs 
into bosons [83"| l84"] and the truncation to only monopole and quadrupole bosons, with the associated 
renormalization of the effective interaction in the truncated space [SoT - 157] . From this perspective, to 
confine to two-body interactions in the boson space, is only a convenient lowest-order approximation. 

The advantages of using higher-order terms with PDS are twofold. First, the algorithms for realizing 
such symmetry structures provide a systematic procedure for identifying and selecting interactions of a 
given order. Having at hand a selection criteria is highly desirable, since, if higher-order terms or new 
degrees of freedom are added, one is immediately faced with the problem of many possible interactions 
and a proliferation of free parameters. Second, Hamiltonians with PDS break the dynamical symmetry 
but retain selected subsets of solvable eigenstates with good symmetry. Such qualities are a virtue, 
since interactions with a PDS can be introduced, in a controlled manner, without destroying results 
previously obtained with a dynamical symmetry for a segment of the spectrum. 

In general, the existence of quantum Hamiltonians with PDS is closely related to the order of the 
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Table 9: Normalized three-boson U(5) tensors. 



n n d r n A t B\ 



n]{n d )(r)n A em 



30000 yl^) 3 

3 110 2 y|(st) 2 4 



3 2 J|st(rftrft)(°) 
3 2 2 2 x /fst(dtdt)(2) 



3 2 2 4 ^|st(rftrft)W 

3 3 1 2 ./A(( rf t d t)(o) d t)W 



3 3 3 1 v / |(( C /tdt)(2) rf t)W 
3 3 3 3 ,/J((dtdt)(2) d t)(^) 



30 1 

3 3 3 4 v /5((rftrft)(2) rf t)W 
3 3 3 6 ,/l((dtdt)(4) rf t)(6) 



interaction among the constituents. IBM Hamiltonians with higher-order terms, exhibiting PDS of 
type II were already encountered in Subsection 13.21 In what follows, we present examples of three-body 
IBM Hamiltonians with U(5) and 0(6) PDS of type I. Work on three-body IBM Hamiltonians with 
SU(3)-PDS of type I, is currently in progress and will be reported elsewhere [55] . 



8.1 U(5) PDS (type I) with three-body terms 

The U(5) dynamical symmetry (DS) chain, U(6) D U(5) D 0(5) D 0(3), and its related basis states, 
\[N](n d )(T)n&LM), were discussed in Section [2TT1 In this case, new terms show up in the DS Hamil- 
tonian at the level of three-body interactions. The DS Hamiltonian and related spectrum now read 



H D s = t 1 n d + t 2 n 2 d + t 3 n 3 d + UC ^ + t 5 n d C (5) +t 6 C 0{3 ) + t 7 n d C ^ , (131a) 
E DS = tin d + t 2 n d + t 3 n d + Ur(r + 3)+t 5 n d T(r + 3)+t (i L(L + l) 

+t 7 n d L(L + l). (131b) 

Terms of the form NCq, with Cq a quadratic Casimir operator of G = U(5), 0(5), 0(3), are included 
in Hds by allowing the parameters £j to depend on N. 

The construction of U(5)-PDS Hamiltonians with three-body terms follows the general algorithm 
by considering operators which annihilate, for example, the U(5) ground state, |[JV],nd = r = L = 0). 
This can be accomplished by means of those U(5) tensors in Table [9] with n d ^ 0. Several families 
of U(5)-PDS Hamiltonians can be defined by identifying specific three-body terms which annihilate 
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additional U(5) basis states. One such family involves the interaction 



V Q = r GlGo + e (gJ# + kJGq) (132a) 
V Q \ [N], n d = r, r, n A = 0, LM) = L = t,t + 1, . . . , 2r - 2, 2r (132b) 

where G^ = [(dW^rff]? 5 with (p,L) = (2,0) (0,2 ), (2,3), (2,4), (4,6) and = S * with 
L = 0, 2, 4. As shown in [89], the states of Eq. ffT32bj) may be projected from states created by acting 
on the vacuum by a product of r operators S m . with > 1 for 2 = 1, . . . , r. Hence, such states are 
guaranteed to contain no three <i-boson states coupled to L = and, therefore, are annihilated by Go- 
The same set of states are annihilated also by K , since all states with = r vanish under the action 
of d ■ d [14J. The remaining eigenstates of V (11321) are mixed with respect to both U(5) and 0(5). 
Clearly, H DS + V exhibits a U(5)-PDS of type I. 

A second family of PDS Hamiltonians involves the interaction 

V 2 = a 2 11^ . C/( 2 ) + J2 e L ^G[-K L + H.c.^j (133a) 

L=0,2,4 

V 2 \[N],n d = r = L = 3) = , (133b) 



where Gl,^ = (— l^Gi-^ and ^ = (—l^Ki-^. The relation in Eq. (1133bj) follows from arguments 



similar to those given after Eq. f )20|) . Other eigenstates of V 2 (I133[) are mixed in the U(5) basis. Clearly, 
H DS + % exhibits a U(5)-PDS of type I. 

A third family of PDS Hamiltonians involves the interaction 

V 3 = r 3 Gl-G 3 + r GlG . (134) 

Both terms in Eq. (1134j) conserve the U(5) quantum number rid, but are not 0(5) scalars. They can 
induce 0(5) mixing subject to At = 2,4,6, and their multipole form involves U(5) generators, some of 
which are not contained in the 0(5) subalgebra. As such, and in accord with the discussion at the end 
of Subsection 13. 2\ H DS + V 3 exhibits U(5)-PDS of type II. Since both terms in V 3 are rotational-scalars 
and are diagonal in n d , it follows that in a given multiplet, those L-states which have a unique 
r-assignment, remain pure with respect to 0(5), and hence are good U(5) eigenstates. For example, 
the states with rid = t, «a = 0, L = 2nd, 2nd — 2, 2nd — 3, 2nd — 5, or states with L = and nd < 5, 
or states with L = 3 and n d < 8, are all eigenstates of V 3 , diagonal in the U(5) basis. In this sense, 
H D s + V 3 exhibits also 0(5)-PDS of type I. 



8.2 0(6) PDS (type I) with three-body terms 

The 0(6) dynamical symmetry (DS) chain, U(6) D 0(6) D 0(5) D 0(3), and its related basis states, 
\\N]{T){r)n^LM) , were discussed in Subsection 12.31 The DS Hamiltonian is given in Eq. (|52p and no 
new terms are added to it at the level of three-body interactions. 

According to the general algorithm, the construction of interactions with 0(6)-PDS of type I requires 
n-boson creation and annihilation operators with definite tensor character in the 0(6) basis: 

B[ n ](a){T)n A em' A™ 5 ] H (r)n A em = (~ 1) (j*[n](ar)(T)n±£,-m) " (135) 

Of particular interest are tensor operators with a < n. They have the property 

B [n s ]{aKr)nA£m \[N}(N)(r)n A LM) =0, a < n, (136) 
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Table 10: Normalized three-boson 0(6) tensors. 



n a t wa t B\ 



n](a)(T)n^£m 



3 3 J^ S t(dtdt)(°) + J^) 3 

3 3 1 2 ^(( t zW)(o) t jt)£> + JlL^ydl 

3 3 2 2 J\s^{SS)^ 

3 3 2 4 W / |st(dtrft)W 

3 3 3 1 y / |((rftrft)(2) rf t)(0) 

3 3 3 3 y / J((rftrft)(2) rf t)© 

3 3 3 4 v /J((dtdt)(2) d t)W 

3 3 3 6 ./^((dW)^)™ 



6* 
7_ 

30 1 
22 ' 

*! 

_5_ 

I((dtrft)(o)dt)W - v 



3 1 ./A s t(dtdt)(°)_ /i( s t)3 



3 110 2 x M((dtdt)(o)rft)W - / i ( a t)2rft 



for all possible values of r,n^,L contained in the 0(6) irrep (N). This is so because the action of 
B\ n &]{c){r)n^im leads to an (N — n)-boson state that contains the 0(6) irreps (£) = (N — n — 2i), % = 
0,1,... which cannot be coupled with (a) to yield (£) = (N), since <r < n. Number-conserving 
normal-ordered interactions that are constructed out of such tensors with a < n (and their Hermitian 
conjugates) thus have \ [N]{N)(r)nALM) as eigenstates with zero eigenvalue (Mj. 

As shown in Subsection I2.3[ there is one two-boson operator Pj = & • & — (s^) 2 , Eq. (|54"|) . with 
a < n = 2, which gives rise to an 0(6)-invariant interaction, P Po, related to the completely solvable 
Casimir operator of 0(6), Eq. fl59|) . On the other hand, from Table [TOj one recognizes two three-boson 
0(6) tensors with a < n = 3 

^[3](l)(l)Q;2m = i^Q^m' ^3](1)(0)0;00 = I^0 S ^ (137) 

and from these one can construct the interactions with an 0(6) PDS. The only three-body interactions 
that are partially solvable in 0(6) are thus P^n^Po and P^n^Po- Since the combination P^hs + fi^Po — 
(N — 2)PjPo is completely solvable in 0(6), there is only one genuine partially solvable three-body 
interaction which can be chosen as P(jn s Po, with tensorial components a = 0, 2. The 0(6)-DS spectrum 



P DS = 4h (N -v + 2)v + Br(r + 3) + C L(L + 1) , (138) 

resembles that of a 7- unstable deformed rotovibrator, where states are arranged in bands with 0(6) 
quantum number £ = iV — 2v, (v = 0, 1, 2, . . .). The 0(5) and 0(3) terms in the dynamical symmetry 
Hamiltonian, H^s (!52j) . govern the in-band rotational splitting. A comparison with the experimental 
spectrum and E2 rates of 196 Pt is shown in Fig. [15] and Table [TT1 The 0(6)-DS limit is seen to provide 
a good description for properties of states in the ground band (S = N). This observation was the 
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Figure 15: Observed spectrum of 196 Pt compared with the calculated spectra of i^DS ([52}) . with 0(6) dynamical 
symmetry (DS), and of -ffpDS (|139j) with 0(6) partial dynamical symmetry (PDS). The parameters in i?os 
(#pds) are h = 43.6 (30.7), B = 44.0 (44.0), C = 17.9 (17.9), and 77 = (8.7) keV. The boson number is TV = 6 
and S is an 0(6) label. Adapted from [24] , 



basis of the claim [90] that the 0(6)-DS is manifested empirically in 196 Pt. However, the resulting fit 
to energies of excited bands is quite poor. The 0^, O3", and O4 levels of 196 Pt at excitation energies 0, 
1403, 1823 keV, respectively, are identified as the bandhead states of the ground (v = 0), first- (v = 1) 
and second- (v = 2) excited vibrational bands [90]. Their empirical anharmonicity, denned by the ratio 
R = E(v = 2)/E(v = 1) - 2, is found to be R = -0.70. In the 0(6)-DS limit these bandhead states 
have r = L = and S = N, N — 2, N — A, respectively. The anharmonicity R = —2/(N + 1), as 
calculated from Eq. fl 1 3 8 j) . is fixed by N. For N = 6, which is the appropriate boson number for 196 Pt, 
the 0(6)-DS value is R = —0.29, which is in marked disagreement with the empirical value. A detailed 
study of double-phonon excitations within the IBM, has concluded that large anharmonicities can be 
incorporated only by the inclusion of at least cubic terms in the Hamiltonian [91] . In the IBM there are 
17 possible three-body interactions [8]. One is thus confronted with the need to select suitable higher- 
order terms that can break the DS in excited bands but preserve it in the ground band. On the basis 
of the preceding discussion this can be accomplished by the following Hamiltonian with 0(6)-PDS [23] 

Hpds = H DS + r, P&sPo, (139) 

where the terms are defined in Eqs. (1321) and (11371) . The spectrum of i^pos is shown in Fig. [15j The 
states belonging to the £ = = 6 multiplet remain solvable with energies given by the same DS 
expression, Eq. (11381) . States with S < 6 are generally admixed but agree better with the data than 
in the DS calculation. For example, the bandhead states of the first- (second-) excited bands have the 
0(6) decomposition E = 4: 76.5% (19.6%), S = 2: 16.1% (18.4%), and S = 0: 7.4% (62.0%). Thus, 
although the ground band is pure, the excited bands exhibit strong 0(6) breaking. The calculated 
0(6)-PDS anharmonicity for these bands is R = —0.63, much closer to the empirical value, R = —0.70. 
It should be emphasized that not only the energies but also the wave functions of the S = N states 
remain unchanged when the Hamiltonian is generalized from DS to PDS. Consequently, the E2 rates 
for transitions among this class of states are the same in the DS and PDS calculations. Thus, the 
additional three-body term in the Hamiltonian f ll39f) . does not spoil the good 0(6)-DS description for 
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Table 11: Observed (EXP) and calculated B(E2) values (in e 2 b 2 ) for 196 Pt. For both the exact (DS) and 
partial (PDS) 0(6) dynamical symmetry calculations, the E2 operator is that of Eq. ([3]) with e# = 0.151 eb 
and x = 0.29. Only the state 0^~ has a mixed 0(6) character. Adapted from |24| . 
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this segment of the spectrum. This is evident in Table [TT] where most of the E2 data concern transitions 
between X = iV = 6 states. Only transitions involving states from excited bands (e.g., the O3" state in 
Table fTTi) can distinguish between Db and PDb. Unfortunately, such interband E2 rates are presently 
poorly known experimentally. Their measurement is highly desirable for further testing the 0(6)-PDb 
wave functions. 



9 PDS and Coupled Systems 

bo far the notion of partial dynamical symmetries was presented in the framework of algebraic models 
involving only one species of constituent particle. It is of great interest to extend this notion to the case 
of coupled systems involving two (or more) species of particles. In this case, the appropriate spectrum 
generating algebra, Gi x G 2 , contains the direct product of the two algebraic structures for systems 1 
and 2. 

An example of such a coupled system is the proton-neutron version of the interacting boson model 
(IBM-2) [8j|92j[93]. The building blocks of the model are monopole and quadrupole bosons, {sL <iL}, 
of proton type (p = tt) and of neutron type (p = v), representing pairs of identical valence nucle- 
ons. Number conserving bilinear combinations of operators in each set comprise the U p (6) algebra 
as in the IBM-1, Eq. (JT]), and bosons of different types commute, bince the separate proton- and 
neutron- boson numbers, N n and N u , are conserved, the appropriate spectrum generating algebra of 
the model is U^(6) x U„(6). bubalgebras can be constructed with the aid of the individual subalgebras, 
U p (5), bU p (3), O p (6), O p (5), O p (3). For instance, for a given algebra G p , with generators g p , there is 
a combined algebra G n + U , with generators g n + g v . The dynamical symmetries of the IBM-2 are ob- 
tained by identifying the lattices of embedded algebras starting with U,r(6) x U„(6) and ending with 
the symmetry algebra O n+u (3). 

A new aspect in coupled systems is the occurrence of states which are not symmetric with respect 
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to interchange of the two constituents. This is clearly seen in the reduction 

U ff (6)xU„(6) D U ff+I/ (6) 

| | . (140) 

[N*]x[N v ] [N U N 2 } 

For a given irrep of U,r(6) x U w (6), characterized by and A^, the allowed irreps of 11^+^(6) are 
[N u N 2 ] = [N n + N u -k, k], where k = 0, 1, . . . , min^, iV^}. States in the irreps with A" 2 7^ (k ^ 0) 
are not symmetric with respect to it and v bosons. One of the successes of the IBM-2 has been the 
empirical discovery of such low-lying mixed symmetry states in nuclei, in which valence protons and 
neutrons move out of phase. A complete listing of all possible partial dynamical symmetries (PDS) of 
the IBM-2 is outside the scope of the present review. In what follows, we present a sample of such 
symmetry structures, illuminating new features of PDS in coupled systems. 

Coupled algebraic structure, G\ x G 2 , can involve also fermionic algebras, as well as Bose-Fermi 
algebras. An example of the latter is the interacting boson-fermion model (IBFM) [5J, used for describ- 
ing odd-mass nuclei and broken fermion-pairs in even-even nuclei. The model incorporates collective 
(bosonic) and quasi-particle (fermionic) degrees of freedom, and the associated spectrum generating 
algebra is Ub(6) x XJ F (m). Here Ub(6) and Up(m) are the boson and fermion algebras respectively, and 
m = Xli(2ji + 1) is the dimension of the single-particle space (ji are the angular momenta of the occupied 
shell-model orbits). Bose-Fermi symmetries correspond to dynamical symmetries of Ub(6) x U F (m). 
Supersymmetry corresponds to a further embedding of the Bose-Fermi symmetry into a graded Lie 
algebra U(6/m) D Ub(6) x Uf(?7j). Partial Bose-Fermi symmetries and partial supersymmetries have 
not been considered in detail so far. There are initial hints that such a structure can occur in the 
IBFM however, an in-depth systematic study is called for. 

9.1 F-spin and selected PDS in the IBM-2 

The proton-neutron degrees of freedom are naturally reflected in the IBM-2 via an SUp(2) F-spin 
algebra [22] with generators 

F + = 4s v + 4 • d u , F_ = (F + )t , F = (N n - N u )/2 . (141) 

These generators commute with the total boson number operator, A" = N n + N u , which is a Ujv(l) 
generator. The basic F-spin doublets are (s^., s\,), and (eft, (ft), with F-spin projection +1/2 (—1/2) 
for proton (neutron) bosons. The algebras SUf(2) x Ujy(l) (I14ip and U n+U (6) f)140p commute and obey 
a duality relationship, in the sense that their irreps are related by F = (Ni — N 2 )/2 = (N n + N u )/2 — k 
and N = Ni + N 2 = N n + N u . In a given nucleus, with fixed N n , N„, all states have the same value 
of F = (N w — N u )/2, while the allowed values of the F-spin quantum number F range from \F \ to 
F max = (N n + N v )/2 = N/2 in unit steps. F-spin characterizes the ir-v symmetry properties of IBM-2 
states. States with maximal F-spin, F = F max , are fully symmetric and correspond to the IBM-1 states 
with only one type of bosons j8]. There are several arguments, e.g., the empirical success of IBM-1, 
the identification of F-spin multiplets [951498] (series of nuclei with constant F and varying Fq with 
nearly constant excitation energies), and weakness of Ml transitions, which lead to the belief that low 
lying collective states have predominantly F = F max States with F < F max , correspond to 'mixed- 
symmetry' states |100j . most notably, the orbital magnetic dipole scissors mode |101] has by now been 
established experimentally as a general phenomena in deformed even-even nuclei |102j . 

Various procedures have been proposed to estimate the F-spin purity of low lying states [99]. In the 
majority of analyses, based on Ml transitions (which should vanish between pure F = F max states), 
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Table 12: Energies (in MeV) of 2 + levels of the ground (g), 7 and (3 bands in F-spin multiplets. The mass 
numbers are A = 132 + 4F. Adapted from |103| . 
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magnetic moments and energy systematics of mixed-symmetry states, the F-spin admixtures in low 
lying states are found to be of a few percents (< 10%), typically 2% — 4% [99J. In spite of its appeal, 
however, F-spin cannot be an exact symmetry of the Hamiltonian. The assumption of F-spin scalar 
Hamiltonians is at variance with the microscopic interpretation of the IBM-2, which necessitates different 
effective interactions between like and unlike nucleons |T7] . Furthermore, if F-spin was a symmetry of 
the Hamiltonian, then all states would have good F-spin and would be arranged in F-spin multiplets. 
Experimentally this is not case. As noted in an analysis [96l|97] of rare earth nuclei, the ground bands 
are in F-spin multiplets, whereas the vibrational /3 bands and some 7 bands do not form good F-spin 
multiplets. The empirical situation in the deformed Dy-Os region is portrayed in Table [T21 and Fig. [TBI 
From Table [121 it is seen that, for F > 13/2, the energies of the L = 2 + members of the 7 bands vary 
fast in the multiplet and not always monotonically. The variation in the energies of the (3 bands is 
large and irregular. Thus both microscopic and empirical arguments rule out F-spin invariance of the 
Hamiltonian. F-spin can at best be an approximate quantum number which is good only for a selected 
set of states while other states are mixed. We are thus confronted with a situation of having 'special 



51 



UJ 



0.4 - 



0.0 



Dy 



Er 



°Hf 



6 0s 



Figure 16: Experimental levels of the ground 7 and j3 bands in an F-spin multiplet F = 6 of rare earth nuclei. 
Levels shown are up to L = 8+ for the ground band, L = 2+ 3+ for the 7 band (diamonds connected by dashed 
lines) and L = 0t,2t for the /3 band (squares connected by dotted lines). Adapted from [103 . 



states' endowed with a good symmetry which does not arise from invariance of the Hamiltonian. These 
are precisely the characteristics of a partial symmetry for which a non-scalar Hamiltonian produces a 
subset of special (at times solvable) states with good symmetry. In what follows we present IBM-2 
Hamiltonians with F-spin as a partial symmetry |1U3] . The construction process is similar to that 
employed in Section for obtaining partially-solvable IBM-1 Hamiltonians. Predictions of F-spin PDS 
are then confronted with empirical data. 

The ground band in the IBM-2 is represented by an intrinsic state which is a product of a proton 
condensate and a rotated neutron condensate with N n and N v bosons, respectively |104j . It depends 
on the quadrupole deformations, (3 P , j p , (p = n, v) of the proton-neutron equilibrium shapes and on the 
relative orientation angles Q between them. For (3 p > 0, the intrinsic state is deformed and members 
of the rotational ground-state band are obtained from it by projection. It has been shown in [105J 
that the intrinsic state will have a well denned F-spin, F = F max , when the proton- neutron shapes are 
aligned and with equal deformations. The conditions (j3 n = /3 v ,j n = j v , Q = 0) are weaker than the 
conditions for F-spin invariance, which makes it possible for a non-F-scalar IBM-2 Hamiltonian to have 
an equilibrium intrinsic state with pure F-spin. Focusing on the most likely situation, namely, aligned 
axially symmetric (prolate) deformed shapes (/3 P — /3, j p — Q — 0), the equilibrium deformed intrinsic 
state for the ground band with F = F max has the form 

\c;K = 0) = \N„N V ) = (NJNjy 1 / 2 ^)^ (&t„)""|0) , 

b{ p = (l + /3 2 r 1/2 (4 + Ko) ; (142) 

where K denotes the angular momentum projection on the symmetry axis. 

The construction of partially-solvable IBM-2 Hamiltonian with F-spin partial symmetry, can be 
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accomplished by means of the following boson-pair operators |103] 



K,o = 4 • 4 - ^ 2 (4) 2 > 4*0,0 = v/2 ( 4 • 4 - /? 2 44 ) 

r{ 2 = V2/3 44 + v / 7(44) (2) , i?{ w))2 = 4 44 + 44 ) + 4H(44) (2) (143) 

wl = (44) (i) (l = i,3), w} = 44 - 44 

The i?^ L pairs (p = 7T, z/) are the same L = 0, 2 pairs of Eq. fl82|) and the 7r-z/ pair, -RLj/) L' com pI e t es 

the set to form an .F-spin vector. Altogether, the R\ L (L = 0,2) are boson pairs with F — 1 and 

(Fq = 1, 0, —1) <->■ [i = 7r, (7rv), i/]. The W 7 ! (L = 1, 2, 3) are F-spin scalar (F = 0) 7r-z/ boson pairs. All 
these operators satisfy 

R iiL ,^\c;K = 0) = 0, (144a) 

W L ^\c;K = 0) = 0, (144b) 

or equivalently, 

Ri,L'/j\c; [N % ], [N v ], F = F max , LM) = 0, (145a) 

W Vil \c- [N v ], [N v ], F = F max ,LM) = . (145b) 

The states in Eq. ( 11451) are those obtained by 7r+l/ (3) projection from the intrinsic state ( 11421) . Since 
the angular momentum projection operator is an F-spin scalar, the projected states of good L will 
also have good F = F max . Following the general algorithm, an IBM-2 Hamiltonian with partial F-spin 
symmetry can be constructed as 

^ = EE A l R U ■ hL + E B ^ W l -W L + C 2 [it^ ))2 • W 2 + H.c] , (146) 

i L=0,2 L=l,2,3 

where Ri,l,h = (— ^YRi,L = The above Hamiltonian is an F-spin scalar only when 

A { [ } = A { 1 ] = A ( l v) (L = 0, 2) and C 2 = 0. Nevertheless, the relations of Eqs. (fllljl - flllg]) ensure that 
it has a solvable zero-energy band with good F-spin for any choice of parameters A®, B^, C 2 and any 
N w , N v . When A®, B L , A { ™ ] B 2 - (C 2 ) 2 > 0, H (USD becomes positive-definite and the solvable states 
form its ground band. We thus have a non-F-spin scalar Hamiltonian with a solvable (degenerate) 
ground band with F = F max . The degeneracy can be lifted by adding to the Hamiltonian 7r+Jy (3) 
rotation terms which produce L{L + 1) type splitting but do not affect the wave functions. States in 
other bands can be mixed with respect to F-spin, hence the F-spin symmetry of H is partial. H trivially 
commutes with F but not with F±. However, [H, F± ]|c; K — 0) =0 does hold and, therefore, H will 
yield F-spin multiplets for members of ground bands. On the other hand, states in other bands can have 
F-spin admixtures and are not compelled to form F-spin multiplets. These features which arise from the 
partial F-spin symmetry of the Hamiltonian are in line with the empirical situation as discussed above 
and as depicted in Table [12] and Fig. [16] It should be noted that the partial F-spin symmetry of H 
holds for any choice of parameters in Eq. (j!46p . In particular, one can incorporate realistic shell-model 
based constraints, by choosing the A 2 (p = 7r, v) terms (representing seniority-changing interactions 
between like nucleons), to be small. For the special choice A^ 2 = C 2 = and B\ = B 3l H of Eq. ( I146j) 
becomes 7r+I/ (5)-scalar which commutes, therefore, with the 0^+^(5) projection operator and hence 
produces F-spin multiplets with good 0^+^(5) symmetry. Such multiplets were reported in the Yb-Os 
region of 7- soft nuclei [98J. 
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The same conditions (f3 p = (3, 7 P = Q = 0) which resulted in F = F max for the condensate of 
Eq. (11 42 p . ensure also F = F max — 1 for the intrinsic state representing the scissors band 

\ SC ;K = 1) = rtjiV.- 1,^-1} , 

rL = &U<i-<i & t- ( 147 ) 

Here Tj c is a F = deformed boson pair whose action on the condensate with (N — 2) bosons produces 
the scissors mode excitation. The states of good L projected from this intrinsic state, \sc; [N n ], [N u ], F = 
Fmax — 1,LM) retain the F-spin quantum number, F = F max — 1. Futhermore, it can be shown that 
the operators R\ Lfi of Eq. (I143j) satisfy 

R i:Ull \sc;K = 1} = , (148a) 
Ri, L ^\sc; [N n ], [N v ], F = F max — 1, LM) = . (148b) 

Consequently, the scissors intrinsic state (11471) and corresponding //-projected states are exact eigen- 
states of the following Hamiltonian, obtained from Eq. (I146P for the special choice C<i = and 
B\ = Z?3 = 2B2 = IB 

H' = S A l R U ■ Ri,L + BMnu ■ (149) 

i L=0,2 

The last term in Eq. (11491) is the Majorana operator [8], with eigenvalues k(N — k + 1) for states with 
F = F max — k. It is related to the total F-spin operator and the quadratic Casimir operator of U vr+Iy (6) 
by M* v = [N(N + 2)/4 — F 2 } = [N(N + 5) - C 2 \U V+U (6)] ]/2. Adding an 0^(3) rotation term we 
obtain the following Hamiltonian with F-spin PDS 

H PDS = H' + C 4(0^(3)] = H DS + Yl a{ l R {l ■ Ri,L ■ (150) 

i L=0,2 

Here Hds contains the Majorana and 0^+^(3) terms, associated with the dynamical symmetry chain 
Utt+j^G) D 7r+I/ (3). Hpds (H50I) has subsets of solvable states which form the K — ground band with 
F = F 

± ± max-, 



\c- [N n ],[N u ],F = F max ,LM) L = 0,2,4,..., 2N (151a) 

E g {L) = CL(L + 1) , (151b) 

and the K = 1 scissors band with F = F max — 1 

\sc; [N K ], [N v ], F = F max - 1, LM) L = 1,2,3, ... , 2N — 1 (152a) 

E SC (L) = BN + CL(L + 1) . (152b) 



It follows that for such Hamiltonians, both the ground and scissors band have good F-spin and have 
the same moment of inertia. The latter derived property is in agreement with the conclusions of a 
comprehensive analysis of the scissors mode in heavy even-even nuclei |106j . which concluded that, 
within the experimental precisions (~ 10%), the moment of inertia of the scissors mode are the same as 
that of the ground band. It is the partial F-spin symmetry of the Hamiltonian (11501) which is responsible 
for the common signatures of collectivity in these two bands. 

The Hamiltonian H' of Eq. (I149p is not F-spin invariant, however, the following relations are satisfied, 
[H' , F] |c; K = 0) = [H', F] \sc; K = 1) =0. This implies that members of both the ground and scissors 
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Table 13: The ratio R = J2 B ( M1 ) t I{Cf,f ) 2 for members of F-spin multiplets. Here ^B(Ml) f denotes 
the experimental summed Ml strength to the scissors mode [1U81I1U9] and Cf,f = (F, Fq;1,0\F — 1,Fq). 
Adapted from [103j . 



Nucleus 


F 


Fo 


££(Ml)t[/4] 


(Cf,f ) 2 


R 


148 Nd 


4 


1 


0.78 (0.07) 


5/12 


1.87 


(0.17) 


148 Sm 




2 


0.43 (0.12) 


1/3 


1 on 

1.29 


(0.36) 


150 Nd 


9/2 


1/2 


1.61 (0.09) 


4/9 


3.62 


(0.20) 


150 Sm 




3/2 


0.92 (0.06) 


2/5 


i on 

2.30 


(0.15) 


154g m 


11/2 


1/2 


2.18 (0.12) 


5/11 


4.80 


(0.26) 


154 Gd 




3/2 


2.60 (0.50) 


14/33 


6.13 


(1.18) 


160 Gd 


7 





2.97 (0.12) 


7/15 


6.36 


(0.26) 


160 Dy 




1 


2.42 (0.18) 


16/35 


5.29 


(0.39) 


162 Dy 


15/2 


1/2 


2.49 (0.13) 


7/15 


5 34 


fO 981 


166 g r 




-1/2 


2.67 (0.19) 


7/15 


5.72 


(0.41) 


164 Dy 


8 





3.18 (0.15) 


8/17 


6.76 


(0.32) 


168 Er 




-1 


3.30 (0.12) 


63/136 


7.12 


(0.26) 


i72 Yb 




-2 


1.94 (0.22) a ) 


15/34 


4.40 


(0.50) 


170 Er 


17/2 


-3/2 


2.63 (0.16) 


70/153 


5.75 


(0.35) 


174 Yb 




-5/2 


2.70 (0.31) 


66/153 


6.26 


(0.72) 



a ) The low value of ^-B(Ml) f for 172 Yb has been attributed to experimental deficiencies |1U2| . 



bands are expected to form F-spin multiplets. For ground bands such structures have been empirically 
established [951498] . The prediction for F-spin multiplets of scissors states requires further elaboration. 
Although the mean energy of the scissors mode is at about 3 MeV [107] . the observed fragmentation 
of the Ml strength among several 1 + states prohibits, unlike ground bands, the use of nearly constant 
excitation energies as a criteria to identify F-spin multiplets of scissors states. Instead, a more sensitive 
test of this suggestion comes from the summed ground to scissors B(M1) strength. The IBM-2 Ml 
operator (L^ — L v ) is an F-spin vector (F — 1,F — 0). Its matrix element between the ground state 
[L = 0g, (F = F max ,F )} and scissors state [L = 1+, (F' = F — l,F )} is proportional to an F-spin 
Clebsch Gordan coefficient Cf,f — {F, F ; 1, 0\F — 1, F ) times a reduced matrix element. It follows that 
the ratio B(M1; 0+ — > l^ c )/ (Cf,f ) 2 does not depend on F and should be a constant in a given F-spin 
multiplet. In Table [TBI we list all F-spin partners for which the summed B(M1) strength to the scissors 
mode has been measured [108. 1109] . It is seen that within the experimental errors, the above ratio is 
fairly constant. The most noticeable discrepancy for 172 Yb (F=8), arises from its measured low value of 
summed B(M1) strength. The latter should be regarded as a lower limit due to experimental deficiencies 
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(large background and strong fragmentation [102] ). These observations strengthen the contention of 
high F-spin purity and formation of F-spin multiplets of scissors states. 

The Hamiltonian H' (11491) depends on (3 through the operators R\ l (I143p . It exhibits additional 
partial symmetries for specific choices of the deformation and/or parameters. Specifically, for f3 = 
\/2, H' has both F-spin and SU 7r+I/ (3) PDS of type I. In such circumstances, the ground (K=0) and 
scissors (K=l) bands have good F-spin and SXJ n+u (3) symmetries: [(A, = [(2N, 0), F max ] and 
[(2N — 2, 1), F = F max — 1], respectively. If in addition, A 2 = A% = A 2 , then also the symmetric-7 
(K = 2), and antisymmetric-7 (K = 2) bands are solvable and have good SU(3) and F-spin symmetries: 
[{2N — 4, 2),F = F max ] and [(2N — 4, 2),F = F max — 1], respectively. In this case, also states of 
the 7-bands will be arranged in F-spin multiplets. At the same time, since the Hamiltonian is not 
F-spin scalar, the ft bands can have F-spin admixtures and need not form F-spin multiplets. As 
noted in [SHUH] an d shown in Table H2] and Fig. [THl such a behaviour is observed for nuclei with 
F — 6, 13/2. For — 1, the ground (K=0) and scissors (K=l) bands have good F-spin and T+I/ (6) 
symmetries: [(<7i, (72}, F] = [(N, 0), F max ] and [(N — l,l),F max — 1], respectively, but the projected 
states are mixed with respect to 0^+^(5) . Consequently, in this case, H' f!149j) has Cv +i ,(6) PDS of 
type III. For A% — -*4. 2 = A ™ = 0, H' is 7r+z/ (5)-invariant. It contains a mixture of terms from 
several chains: U 7r (5) x U u (5), 0^(6) x Cv(6) and XJ W+U (6), all sharing a common 7r+ly (5) D 7r+Jy (3) 
segment. In such circumstances, H' has a partially-solvable 0^+^(5) PDS of type II. Such a PDS was 
used in jllUj to obtain an extended Ml sum rule for excited symmetric and mixed-symmetry states, 
and apply it to 94 Mo. 

A new aspect that can occur in an algebraic description of coupled systems, is the situation in which 
the set of Casimir operators in a given chain of subalgebras of G, may not be sufficient to express 
the most general Hamiltonian constructed from the generators of G. Such a scenario was considered 
in [11 lj in connection with the U n+lJ (5) chain of the IBM-2, and shown to be associated with a partial 
dynamical symmetry. The U n+U (5) limit of the IBM-2 corresponds to the chain 

U w (6) x U„(6) D IW6) D IW(5) D O w+I/ (5) D O w+v (3) 

I I ill. (153) 

[N v ]x[N„] [N-k,k] K,n 2 } (n,r 2 ) a L 

The total number of bosons is N = + N v and their F-spin is F = N/2 — k. The states conserve also 
the total number of (i-bosons, rid = rii + "-2, and their separate F-spin, Fa = (rii — n^jl. Apart from 
N p - and ^-dependent terms, the most general one- and two-body Hamiltonian which has a U 7r+iy (5) 
dynamical symmetry (DS) is given by 

Hds = eC 1 [\J 7r+u (5)] + r ] (C 1 [V 7T+ ,(5)}) 2 + AC 2 [V 7r+u (5)} 

+Bd 2 [O v+v (S)] + Cd 2 [O v+v (3)] + aM vt , , (154) 

where C P [G] denoted the p-th order Casimir operator of G. However, it is not the most general Hamil- 
tonian constructed from the generators of U n+U (5). To obtain the latter, another independent operator 
must be added to Hos which is not a Casimir operator of a subalgebra. A simple choice of such an 
operator can be [11 lj 

Vi = (155) 

where W\p is the (F = 0, L = 1) boson-pair defined in Eq. (11431) . The latter transforms as a (ti,Tz) = 
(1,1) tensor under 0^+^(5). Consequently, V\ has components with (ti,t 2 ) = (2,2) © (0,0), hence 
breaks the 7r+l/ (5) symmetry. The 11^+^(6) irrep [N — k, k] with k ^ contains 0^+^(5) irreps (n, r 2 ) 
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with t 2 7^ 0, which can be admixed by this term. Nevertheless, V\ has a subset of zero-energy solvable 
states with good 0^+^(5) symmetry These are the 11^+^(5) basis states with F d = n d /2, which are 
annihilated by 

W^liK] x [N v ]; [N-k,k],{n d ,0},(r,Q),n A ,LM) = 0. (156) 

The interaction V\ can be added to Hds (1154j) to obtain the following Hamiltonian with 0^+^(5) PDS 
of type I 

Hpds = Hds + ZiWI-Wx . (157) 
Hpds breaks the 11^+^(5) DS but retains a subset of solvable U^ +i/ (5) basis states with known eigenvalues 
\[N n ] x [N„\; [N-k,k},{n d ,0},(r,0),n A ,LM) 

Ep DS = en d + V n 2 d + An d (n d + 4) + Bt(t + 3) + CL(L + 1) + ak(N - k + 1) . 1 °° } 

Hpds (H57p exhibits also a U 7r+! ^(6) D XJ n+u (5) PDS of type II, since the remaining eigenstates pre- 
serve the quantum numbers of \J n+u (6), \J n+v (5) and 0^+^(3) but not of 7r+!/ (5) in the chain (I153p . 
V\ (I155P has additional zero-energy eigenstates with F d < n d /2 which break, however, the 0^+^(5) sym- 
metry |lllj . These solvable states lead to additional PDS of the Hamiltonian (I157p . provided B = in 
Eq. f ll54p . In general, PDS associated with vanishing eigenvalues of V\ can explain the simple regulari- 
ties in the spectra of the generalized Majorana operator, observed in an IBM-2 analysis of Pd and Ru 
nuclei [TT2] . 

10 PDS in Fermion Systems 

Partial symmetries are not confined to bosonic systems. The proposed algorithms for constructing 
Hamiltonians with PDS do not rely on the statistics of the constituents, hence can be implemented for 
both bosons and fermions. Identifying partial symmetries in fermion systems can proceed in two ways. 
The first approach relies on a mapping of a bosonic Hamiltonian, which posses a partial symmetry, 
into its fermionic counterpart. If the bosonic generators of the spectrum generating algebra can be 
mapped into fermionic generators of the same algebra, then both Hamiltonians will exhibit the same 
type of partial symmetry. This approach was demonstrated in |113j for schematic U(2) x U(2) Lipkin- 
type models. A second approach relies on a direct construction of fermion Hamiltonians with partial 
symmetries. In what follows, we demonstrate this approach by identifying fermionic PDS related to 
properties of the quadrupole-quadrupole interaction in the framework of the symplectic shell-model [1141 
1115] and to properties of seniority-conserving and non-conserving interactions in a single j shell |116f - 
122J. Such findings constitute a first step towards understanding the microscopic origin of PDS in 
nuclei. 

10.1 PDS in the symplectic shell model 

The symplectic shell model (SSM) [5| I123] is an algebraic, fermionic, shell-model scheme which includes 
multiple 2fajj one-particle one-hole excitations. The scheme is based on the symplectic algebra Sp(6,R) 
whose generators Af^\ C^' and H have good SU(3) [superscript (A,/i)] and 0(3) [subscript 

£, m] tensorial properties. The Af^ = (— ^Y~ m {Af'^ m )^], £ = or 2, create (annihilate) 2fkj 

excitations in the system. The Cj^, £ = 1 or 2, generate a SU(3) subgroup and act only within one 
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[4 to] 



(k+4,0) (X+2,1) (A.,2) (A,- 1,1) (k-2,0) (k-2,3) (k-4,4) 



[2 to] 



{0^2,02,03} {(^+1,02+1,03} {a 1 ,a 2 +2,a 3 } 



[Oto] 



(A.+2,0) (A.,1) (X-2,2) 



{a 1 ,a 2 ,a 3 ) 



{A x A}( 4 °) 
{A x A}(° 2 ) 



A(20) 



Figure 17: Basis construction in the symplectic model. SU(3)-coupled products of the raising operator A- 2 ^ 
with itself act on an Elliott starting state with (A CT , p a ) = (A, 0) ({01, o"2> 03 = 02}) to generate symplectic 2hu, 
4hu}, ... excitations. Also shown are the SU(3) labels (A, p) and quanta distributions {u)\, 002, W3} for some 
excited states. Adapted from [114] . 



harmonic oscillator (h.o.) shell (V^C^ = Qimi the symmetrized quadrupole operator of Elliott, which 
does not couple different h.o. shells [I], and = L m , the orbital angular momentum operator). 

The harmonic oscillator Hamiltonian, H , is a SU(3) scalar and generates U(l) in U(3) = SU(3) x 
U(l). A fermion realization of these generators is given in |124j . The model fully accommodates 

the action of the collective quadrupole operator, Q 2 m = y ^ ^2 S r 2 Y 2m (f s ), which takes the form, 

A basis for the symplectic model is generated by applying symmetrically coupled products of the 2ftw 
raising operator A^ 20 ^ with itself, to the usual Ohio many-particle shell-model states. Each Ofauj starting 
configuration is characterized by the distribution of oscillator quanta into the three cartesian directions, 
{<7i, 02, 03} (01 > 02 > 03), or, equivalently, by its U(l)xSU(3) quantum numbers N a (A CT , p a ). Here 
A CT = 0i — 2 , p a = &2 — 0-3 are the Elliott SU(3) labels, and N a = o\ + 2 + 03 is related to the eigenvalue 
of the oscillator number operator. Each such set of U(3) quantum numbers uniquely determines an 
irrep of the symplectic group, since it characterizes a Sp(6,R) lowest weight state. The product of 
N/2, N = 0,2,4,..., raising operators is multiplicity-free and generates Nhw excitations for 

each starting irrep N a (X a ,p a ). Each such product operator p 7V ( A ™'^") ) labeled according to its SU(3) 
content, (X n ,p n ), is coupled with \N a (X a , p a )) to good SU(3) symmetry p(X,p), with p denoting the 
multiplicity of the coupling (X n , p n ) x (X a ,p a ). The quanta distribution in the resulting state is given 
by {coi, CO2, ^3}, with N a + N = Ui + U2 + coi > 0J2 > ^3, and A = u\ — u>2, p = U2 — CO3. The 
basis state construction is schematically illustrated in Fig. [T7] for a typical Elliott starting state with 
(X a , p a ) = (A,0). To complete the basis state labeling, additional quantum numbers a = kLM are 
required, where k is a multiplicity index, which enumerates multiple occurrences of a particular L value 
in the SU(3) irrep (A, p) from 1 to K™ ax (X, p) = [(A + p + 2 - L)/2] - [(A + 1 - L)/2\ -[(p + 1 - L)/2], 
where [. . .] is the greatest non- negative integer function [125]. The orthonormal SU(3) basis employed 
is that of Vergados |30j, however, for convenience, the running index k — 1,2, ... , K™ ax is used instead 
of the usual Vergados label, \. The dynamical symmetry chain and the associated quantum labels of 
the above scheme are given by |123j : 

Sp(6,R) D U(3) D SO(3) 

I I . (159) 

N a (X a ,p a ) N(\ n ,p n )p N^Xu.fiu) k L 
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Figure 18: Energy spectra for 12 C. Comparison between experimental values (left), results from a symplectic 8Huj 
calculation (center) and a PDS calculation (right). K=0i indicates the ground band in all three parts of the hgure. In 
addition, resonance bands dominated by 2fku excitations (K=2 1; 02, 11,03), 4Huj excitations (K=4i), and 6huj excitations 
(K=6i) are shown for the Sp(6,R) and PDS calculations. Additional mixed resonance bands (not shown), dominated by 
Afjw and Qfvjj excitations, exist for this nucleus. The angular momenta of the positive parity states in the rotational bands 
are L=0,2,4,. . . for K=0 and L=K,K+l,K+2, . . . otherwise. Bands which consist of pure-SU(3) eigenstates of the PDS 
Hamiltonian (|164j) are indicated. Adapted from [115] . 

The following SSM Hamiltonians which has SU(3) partial symmetry have been proposed [114]lll5j 
£(A>,/3 2 ) = P A B + p 2 A 2 -B 2 

= ^ (93*1(3) - 9<7sp( 6 ) + 3# 2 - 3QH ) + (A> - /3 2 )ioA) . (160) 

Here Ag m = AfJ^, Bn m = Bfjf* and the Casimir operators, Cg, conform with the conventions used 
in [114LI115] . For /3q = (3 2 , the Hamiltonian is an SU(3) scalar which is diagonal in the dynamical 
symmetry basis (11591) . For /3q = —5ft 2 , the Hamiltonian transforms as a (2, 2) tensor under SU(3). Thus, 
in general, H((3 , (3 2 ) is not SU(3) invariant, however, it exhibits partial SU(3) symmetry. Specifically, 
among the eigenstates of H((3 ,(3 2 ), there exists a subset of solvable pure-SU(3) states, \4>lm{N)), the 
SU(3)dO(3) classification of which depends on both the Elliott labels (A^/v) of the starting state 
and the symplectic excitation N. In general, it is found that all L-states in the starting configuration 
(N = 0) are solvable with good SU(3) symmetry (X a ,fJ, a ). For excited configurations, with N > (N 
even), one can distinguish between two possible cases: 

(a) A CT > fi a : the pure states belong to (A,//) = (A CT — N,fj, a + N) at the Nfiw level and have 
L = n a + N, ii a + N + 1, . . . , X a - N + 1 with N — 2,4, . . . subject to 2N < (X a - \j, a + 1). 

(b) A CT < the special states belong to (A,/i) = (A CT + N,fj, a ) at the Nfkv level and have L = 
X a + N,X a + N+l,...,X a + N + fx a with N — 2,4, — 

To prove the claim, it is sufficient to show that Bq annihilates the states in question 

%\4> LM {N)) = 0. (161) 

For N = this follows immediately from the fact that the Ohu starting configuration is a Sp(6,R) 
lowest weight which, by definition, is annihilated by the lowering operators of the Sp(6,R) algebra. 
The latter include the generators B^\ For > 0, let {<Ji,<J 2 ,a 3 } be the quanta distribution for 
a Ohw state with X a > Adding quanta to the 2-direction yields a Nhu state with quanta 
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Table 14: B(E2) values (in Weisskopf units) for ground band transitions in 12 C. Compared are several symplec- 
tic calculations, PDS results, and experimental data. Q denotes the static quadrupole moment of the L w = 2f 
state in units of eb. PDS results are rescaled by an effective charge e*=1.33 and the symplectic calculations 
employ bare charges. Adapted from [115J. 



Transition 




Model B(E2) [W.u.] 




B(E2) [W.u.] 


Ji ->■ Jf 


2Huj 




6hu 


8hu 


PDS 


Exp 




4.65 


4.65 


4.65 


4.65 


4.65 


4.65 ± 0.26 


4^2 


4.35 


4.27 


4.24 


4.23 


4.28 


n/a 


Q [eb] 


0.059 


0.060 


0.060 


0.060 


0.058 


0.06± 0.03 



Table 15: B(E2) values (in Weisskopf units) for ground band transitions in 20 Ne. Compared are several 
symplectic calculations, PDS results, and experimental data. The static quadrupole moment of the 2f state is 
given in the last row. PDS results are rescaled by an effective charge e*=1.95 and the symplectic calculations 
employ bare charges. Adapted from [115] . 



Transition 
Ji -> Jf 




Model B(E2) [W.u.] 




B(E2) [W.u.] 
Exp 


2hu 




6hu 


8^; 


PDS 


2^0 


14.0 


18.7 


19.1 


19.3 


20.3 


20.3 ± 1.0 


4^2 


18.4 


24.5 


24.6 


24.5 


25.7 


22.0 ± 2.0 


6^4 


17.1 


22.3 


21.5 


20.9 


21.8 


20.0 ± 3.0 


8^6 


12.4 


15.2 


13.3 


12.4 


12.9 


9.0 ± 1.3 


Q [eb] 


-0.14 


-0.16 


-0.16 


-0.16 


-0.17 


-0.23 ± 0.03 
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Figure 19: Energy spectra for 20 Ne. Experimental ground band (K=0i) energies (left), compared to theoretical 
results for both the ground band and 2fkj resonances (K=02, li, 21,03) for a symplectic 8hw calculation (center) 
and a PDS calculation (right). Rotational bands which consist of pure eigenstates of the PDS Hamiltonian are 
indicated. Adapted from [115J. 



distribution {cr 1 , o 2 + N, cr 3 }, that is, (A, fx) = (X a — N, fM a + N). Acting with the rotational invariant B 
on such a state does not affect the angular momentum, but removes two quanta from the 2-direction, 
giving a (N — 2)fko state with (A',//) = (A CT — iV + 2, fi a -\- N — 2). (The symplectic generator B 
cannot remove quanta from the other two directions of this particular state, since this would yield a 
state belonging to a different symplectic irrep.) Comparing the number of L occurrences in (A,/i) and 
(A',//), one finds that as long as X a - N + 1 > + N, A L (N) = K™ ax (\,/j,) - K^ ax (X'^') = 1 for 
L = ^ + N, fi a + N + 1, . . . , A CT - JV + 1, and Ai(N) = otherwise. Therefore, when A L (N)=1, a 
linear combination |0i(iV)) = c K \Nhuj(X a — N,/j la + N)kLM) exists such that B \(()l(N)) = 0, and 
thus the claim for family (a) holds. The proof for family (b) can be carried out analogously. Here the 
special irrep (A,/i) = (A CT + N,fi a ) is obtained by adding quanta to the 1-direction of the starting 
configuration. In this case there is no restriction on N, hence family (b) is infinite. 

The special states have well denned symmetry with respect to the chain fl 1 5 9 [) and the condition of 
Eq. ( 11611) ensures that they are solvable eigenstates of ^(/^o,^) ( 11601) with eigenvalues ^(A^ = 0) = 
for the Oftw level, and 

N 3 
E(N) = p 2 j{N a -X a + fi a -6+-N) (A CT >/v) (162a) 

N 3 
E(N) = p 2 -(N a + 2X a + fi a -3 + -N) (A CT </v) (162b) 

for N > 0. All Ohu states are unmixed and span the entire (A CT ,/i CT ) irrep. In contrast, for the excited 
levels (N > 0), the pure states span only part of the corresponding SU(3) irreps. There are other states 
at each excited level which do not preserve the SU(3) symmetry and therefore contain a mixture of SU(3) 
irreps. The partial SU(3) symmetry of ^(^0,^2) is converted into SU(3)-PDS of type I by adding to 
it 0(3) rotation terms which lead to L(L+l)-type splitting but do not affect the wave functions. The 
solvable states then form rotational bands and since condition (11611) determines their wave functions, 
one can obtain analytic expressions for the E2 rates between them [115j . 

The Hamiltonian of Eq. (11601) . with SU(3) partial symmetry, has a close relationship to the collective 
quadrupole- quadrupole interaction, 

Q 2 • Q 2 = 3(C 2 + A 2 + B 2 ) ■ {C 2 + A 2 + B 2 ) 

= H(/3q = 12, 02 = 18) + const — 3L 2 + {terms coupling different h.o. shells} . (163) 
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Fi gure 20; Decompositions for calculated 2^ states of 2 "Ne (left) and 12 C (right). Individual contributions from the 
relevant SU(3) irreps at the Ohu and 2hu> levels are shown for both a symplectic 8hoj calculation (denoted Qi ■ Qi) and 
a PDS calculation. In addition, the total strengths contributed by the Nhui excitations for N > 2 are given for the 
symplectic case. Adapted from |115j . 



Here const = 6C'sp(6) — 2Hq + 3AHq is fixed for a given symplectic irrep N a (X a , fj, a ) and Nfrw excitation. 
Although H(/3 , /3 2 ) does not couple different harmonic oscillator shells, it breaks the SU(3)-symmetry 
and, due to the above relation, exhibits in-shell behavior similar to that of Q2 ■ Qi- To study this 
connection further, and to illustrate that the PDS Hamiltonians of Eq. ( 11601) are physically relevant, 
the formalism was applied to oblate prolate and triaxially deformed nuclei, by comparing the properties 
of the following SU(3)-PDS Hamiltonian 

Hpds = h(N) +f#(A> = 12, & = 18) + l2 L 2 + l4 L 4 (164) 

to those of the symplectic Hamiltonian 

^sp( 6) = Ho ~ VQ2 ■ Q2 + d 2 L 2 + d,L 4 . (165) 

Here the function h(N), which contains the harmonic oscillator term Ho, is simply a constant for a 
given Nfyjj excitation. Details of the calculations can be found in [115] . 

The leading Sp(6,R) irrep for the oblate nucleus 12 C is N a (X a ,fi a ) = 24.5(0,4). In this case, 
Hpds ( 11641) has the pure-SU(3) states of family (b) as solvable eigenstates. In particular, all Ohw states 
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Figure 21: Energy spectra for 24 Mg. Energies from a PDS calculation (bottom) are compared to 
symplectic results (top). Both Ohco dominated bands (K=0i,2i,4i) and some 2hco resonance bands 
(K=02, 03, 04, 22, 23, 42, 43, 61) are shown. The K=0i,2i,4i (61) states are pure (approximately pure) in the 
PDS scheme. Experimental energies of the ground and 7 bands are shown on the left. Adapted from [115J. 



are pure (A CT , /v) = (0, 4) states, and at 2hw a rotational band with good SU(3) symmetry (A, fi) = (2, 4) 
and L = 2,3,4,5,6 exists. Similarly, at Afajj there are pure-SU(3) bands with (A,//) = (4,4) and 
L — 4,5, 6, 7, 8, at 6hw with (A, /x) = (6, 4) and L = 6,7, 8, 9, 10, etc. On the other hand, the leading 
Sp(6,R) irrep for the prolate nucleus 20 Ne is N a (X a ,fj, a ) = 48.5(8,0). In this case, the solvable pure- 
SU(3) eigenstates of Hpns are those of family (a). They form a K=0i L = 0,2,4,6,8 rotational band 
with (A,/i) = (8,0) at Ohw, a K=2 t L = 2,3,4,5,6,7 band with (A,/i) = (6,2) at 2kj, and a K=4j 
L = 4, 5 'band' with (A, fx) = (4,4) at Atvjj. There are no other pure PDS states at higher levels of 
excitation. As shown in Fig. [TSIand Table [T41 for 12 C, and in Fig. [T9land Table [TBI for 20 Ne, the resulting 
energies and transition rates in the PDS and symplectic approaches are similar and converge to values 
which agree with the data. The PDS calculations of B(E2) values, however, require an effective charge. 
A better measure for the level of agreement between the PDS and symplectic results is given by a 
comparison of the eigenstates. From Fig. [2D] we observe that, although Hsp{6) (H65p [Hpos (1 164}) ] does 
(does not) mix different major oscillator shells, the Nhoj level to which a particular PDS band belongs 
is also dominant in the corresponding symplectic band. Furthermore, within this dominant excitation, 
eigenstates of Hsp(6) an d H PDS have very similar SU(3) structure, that is, the relative strengths of the 
various SU(3) irreps in the symplectic states are approximately reproduced in the PDS case. This holds 
for the ground and excited bands in both nuclei, with the exception of the K=02 resonance band in 
20 Ne, where significant differences appear in the structure of the wave functions. Thus the SU(3)-PDS 
Hamiltonian captures much of the physics of the Q2 ■ Q2 interaction. 

Calculations were also performed [115] for the triaxially-deformed nucleus, 24 Mg, whose leading 
Sp(6,R) irrep is N a (X a , fi a ) = 62.5(8,4). Since now both A CT 7^ and \i a 7^ 0, the symplectic Hilbert 
space has a very rich structure. In this case, the solvable pure-SU(3) states are those of family (a) and 
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Figure 22: Decompositions for calculated 17 = 6 + states of 24 Mg. Eigenstates resulting from the symplectic 6Hu> 
calculation (denoted Qi -Q-i) are decomposed into their Ohuj, 2huj, 4huj, and 6Huj components. At the Ohio and 2fku levels, 
contributions from the individual SU(3) irreps are shown, for higher excitations (N > 2) only the summed strengths are 
given. Eigenstates of the PDS Hamiltonian belong entirely to one Nhjj level of excitation, here Ohaj or 2tko. Contributions 
from the individual SU(3) irreps at these levels arc shown. Members of the K=0i, 2i, 4i bands are pure in the PDS scheme, 
and K=6i states are very nearly (> 99%) pure. Adapted from |115j . 



include three rotational K=0,2,4, bands at Ohu with (A,/i) = (8,4), and at 2hu a (short) rotational 
K=6 band with L—6,7, which belongs to the irrep (A, fi) = (6, 6). Figure |2~T1 compares the experimental 
spectrum of 24 Mg with energies obtained with 6hw symplectic and PDS calculations. Both the PDS and 
symplectic Hamiltonians were supplemented with small cubic and quartic SU(3)-conserving interactions 
to account for .fT-band splitting. These extra terms break the partial symmetry slightly (< 1%) in the 
K = 6i band, as can be inferred from the eigenstate decompositions plotted in Fig. [22j As in the 
previous examples, the eigenstates of the PDS and symplectic Hamiltonians are found to have very 
similar structures. The structural differences that do exist are reflected in the very sensitive interband 
transition rates [115] . 

The boson Hamiltonian, H(ho, h%) of Eq. (1341) . and the fermion Hamiltonian, H(/3q, 02) of Eq. (I160p . 
have several features in common. Both display partial SU(3) symmetry, they are constructed to be 
rotationally invariant functions of (A, fi) = (2,0) and (A,/i) = (0,2) SU(3) tensor operators, and both 
contain components with (A,/i) = (0,0) and (2,2). Both Hamiltonian have solvable pure-SU(3) eigen- 
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states, which can be organized into rotational bands. The ground bands are pure in both cases, and 
higher-energy pure bands coexist with mixed-symmetry states. There are, however, several significant 
differences between the bosonic and fermionic PDS Hamiltonians. For example, the ground band of 
the bosonic Hamiltonian H(h ,h 2 ), Eq. ( |34|) . is characterized by (A,/i) = (2N, 0), i.e., it describes an 
axially-symmetric prolate nucleus. As mentioned at the end of Subsection 12.21 it is also possible to 
find an IBM Hamiltonian with partial SU(3) symmetry for an oblate nucleus. It can be shown that 
these two cases exhaust all possibilities for partial SU(3) symmetry with a two-body Hamiltonian in 
the IBM with one type of monopole and quadrupole bosons. In contrast, the fermionic Hamiltonian 
H(/3o, 02) (H60p can accommodate ground bands of prolate [(X a , 0)], oblate [(0, p a )], and triaxial [(X a , p a ) 
with Ao- 7^ 0, p a 7^ 0] shapes. Another difference lies in the physical interpretation of the excited solv- 
able bands. While these bands represent j k excitations in the IBM, they correspond to giant monopole 
and quadrupole resonances in the fermion case. Furthermore, whereas the pure eigenstates of the IBM 
Hamiltonian H(ho, h 2 ) can be generated by angular momentum projection from intrinsic states (|37|) . a 
similar straightforward construction process for the special eigenstates of the symplectic Hamiltonian 
H(j3 , /3 2 ) has not been identified yet. The situation seems to be more complicated in the fermion 
case, which is also reflected in the fact that if^o,/^) has two possible families of pure eigenstates, one 
finite, the other infinite. The association of the special states to one or the other family depends on 
the Ohu symplectic starting configuration. Thus, in spite of similar algebraic structures, the bosonic 
and fermionic PDS Hamiltonians involve different mechanisms for generating the partial symmetries in 
question. 

10.2 PDS and seniority 

The notions of pairing and seniority for n identical nucleons occupying a single j shell, are conveniently 
encoded in the quasi-spin formalism [3]. The latter is based on an SUg(2) algebra whose generators are 

S t = \f^ (4 Q ])S Q) > S i = ( S ty > S°=~(n-fi). (166) 

Here Sj~ creates a pair of fermions with angular momentum J = 0, h = J2 m a ]m a jm is the number- 
operator and fl = (2j + l)/2. The basic quasi-spin doublet is (aj m , —dj m ), with dj m = (—iy +m aj- m . 
The seniority quantum number, v, refers to the number of nucleons which are not in zero-coupled J = 
pairs. The SUs(2) quantum numbers, (S, So) are related to n and v by 

S = (n-v)/2 , S = (n-n)/2 . (167) 

The operator S~ annihilates states with (S, S Q = —S) corresponding to states in the j v configuration 
with seniority v and n = v. 

The relevant chain of nested algebras for n fermions in a single j shell is |126[I127] 

U(2j + 1) D Sp(2j + 1) D SUj(2) 

I I 4 (168) 

[l n ] v p J 

where p is a multiplicity index. The generators of the indicated algebras are 

U(2j + 1) = {f. i L) } , Sp(2j + 1) = {f W Lodd} , 

SUj(2) = {J m = cjfV} , Cj = V / j(j + l)(2j + l)/3 , (169) 
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where 



f^=(a)~a^ m L = 0,l,2,...,2j . (170) 

The unitary symplectic algebra, Sp(2j+1), and quasi-spin algebra, SUg(2), commute. The duality 
relationship between their respective irreps, v and S, is given in Eq. (11671) and their quadratic Casimir 
operators are related by C*2[Sp(2j + 1)] = Q(Q + 2) — AS 2 . The dynamical symmetry Hamiltonian 
associated with the chain (J168p is given by 

H DS = pn + an(n-l) + aC 2 [Sp(2j + l)] + bC 2 [SXJj(2)} , (171a) 
Eds = (3n + an(n-l) + av(2tt + 2-v) + bJ(J + l) , (171b) 

and E D s are the corresponding eigenvalues. 

The fi L) operators of Eq. flTTOj) , with L odd, are quasi-spin scalars. Accordingly, the most gen- 
eral number-conserving rotational invariant Hamiltonian with seniority- conserving one- and two-body 
interactions, acting in a single- j shell, can be transcribed in the form [T7| 

H = /3n + an(n- 1) + A L f (L) ■ f (L) (172a) 

Lodd 

E n ,v, P ,j = 13 n + an(n - 1) + Z V)P)J . (172b) 

For the special choice, Xl = 2a + he 1 - 5l,i, the above Hamiltonian reduces to the dynamical symmetry 
Hamiltonian of Eq. (I171al) . In general, eigenstates of H (11721) are basis states, \n,v, p, JM) , of the 
chain (I168P with eigenvalues E n)V)P) j and wave functions of the form 

/ A \ (n— v)/2 

\n,v,p,JM)otiS+) \n = v,j,v,p,JM) . (173) 

Since the last term in Eq. (I172a|) is a quasi-spin scalar, its eigenvalues, Z VtPt j, are independent of n. 
Consequently, for a given n, energy differences E n j tVtPt j — E n j iV / y t j> , are independent of n. Conservation 
of seniority does not, however, imply solvability. In general, eigenstates and eigenvalues of H, Eq. (I172p . 
must be obtained from a numerical calculation. Nevertheless, it has been shown that some multiplicity- 
free eigenstates of H, i.e., with unique n, v, J assignments, are completely solvable and closed algebraic 
expressions for their energies have been derived [116U1 17j . These include states with {y = 2, J), (v, J m ax) 
and (v, J max — 2), where J ma x = v(2j + l — v)/2. As an example, for n even, the ground state of H (I172a[) . 
with v = J = 0, is analytically solvable with energy 

E n ,j,v=o,j=o = /3n + an(n-i) . (174) 

The above expression is identical to the ground-state energy of the dynamical symmetry Hamiltonian, 
Eq. (I17ip . Since all eigenstates carry the seniority quantum number v , the Hamiltonian H of Eq. (11721) 
exhibits PDS of type II with an added feature that some multiplicity-free states are analytically solvable. 

As is well known [17] , seniority remains a good quantum number for any two-body interaction acting 
within a j shell when j < 7/2, but it need not be conserved for j > 9/2. Recently, it has been shown 
that it is possible to construct interactions that in general do not conserve seniority but which have 
some solvable eigenstates with good seniority |118[I120] . Specifically, for j = 9/2, n = 4, there are two 
independent states (p = 1, 2) with seniority v = 4 and J = 4, 6. For each such J value, there exists 
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one particular combination, which is completely solvable with good seniority, v = 4, for any two-body 
interaction in the j = 9/2 shell. The energies of the solvable states are given by [120] 

114 



E[(9/2)\v = 4,J 



£[(9/2)> = 4,J 



68 13 
4] = -u 2 + u A + -u G 



55 



"8 



336 



19 12 

6 = — v 2 H va + vf, , 

j 11 2 ^13 6 143 



^8 



(175a) 



(175b) 



where v\ = ((9/2) 2 ; A|V"|(9/2) 2 ; A) are the matrix elements of an arbitrary rotational-invariant two-body 
interaction, V. The indicated states retain their structure and are completely solvable, independent of 
whether the interaction conserves seniority or not. It follows that the most general one- and two-body 
rotational- invariant Hamiltonian in the j = 9/2 shell, exhibits PDS of type I. The E2 matrix element 
between the two solvable states is interaction- independent, and the corresponding B(E2) value is given 

by m\ 



B(E2; (9/2) 4 , v = 4, J = 6 -> (9/2) 4 , v = 4, J = 4) = 

2D9475 

B{E2- (9/2) 2 , J = 2 -+ (9/2) 2 , J = 0) . 



(176) 



176468 

This again defines a parameter-independent relation between a property of the two- and four-particle 
system. Some properties of these solvable states can be attributed to properties of certain coefficients of 
fractional parentage [119H122] . It has been suggested that this partial seniority conservation may shed 
some new light on the existence of isomers in nuclei with valence neutrons or protons predominantly 
confined to the g 9 / 2 or hg/ 2 shell [120] . 

In most medium-mass and heavy nuclei, the valence nucleons are distributed over several non- 
degenerate j orbits in a major shell. This situation can be treated within the generalized seniority 
scheme [128] - based on more general pair-operators with J = 0, 2 



3,3 



The generalized seniority conditions [128H130] are 

F|0) = 0, 
|0> = V 5t|0> , 

|o> = ^l|o> , 



H, 



H , 5 f 
H, 



Dt 



A &Dt 



where |0) is the doubly-magic core state. These relations imply that the Hamiltonian has 
ground state with J = and generalized seniority v g — 



\n = 2N,v g 



0,J = 0) oc (&) N \0) 



-Eat,« 9 =o,j=o 



NV + -N(N- 1) A 



and a solvable excited state with J = 2 and generalized seniority v„ 



\n = 2N,v g = 2,J = 2) oc (S ] ) N 1 D^) 



En,v b =2,J=2 — -EtV,-!; 9 =0,J=0 + V 2 — Vq 



(177a) 
(177b) 

(178a) 
(178b) 

(178c) 

a solvable 

(179a) 
(179b) 



(180a) 
(180b) 
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From the point of view of symmetry, the monopole pair-operators S^, S, and [S^ , S~\ /2, of Eq. f 1 1 77aj) . 
with unequal ctj, do not form an SU(2) algebra. Nevertheless, an Hamiltonian H satisfying rela- 
tions (I178p has selected solvable states, i.e., it is partially solvable. It should be noted that the energy 
and wave function of the v g = ground state (I179P have the same form as in the single- j dynamical 
symmetry expressions, Eqs. (I173l) - fll74p and, by construction, both spectra exhibit linear two-nucleon 
separation energies and constant 2 + -0 + spacings. Furthermore, the Sp(2j+1) generators (I169p . Tm 
with L odd and any j, annihilate the state of Eq. (I179p . Therefore, the v g = ground state is invariant 
under JT. Cg>Sp(2j + 1), even though the latter is not a symmetry of the Hamiltonian. 



11 Concluding Remarks 

The notion of partial dynamical symmetry (PDS) extends and complements the fundamental concepts 
of exact and dynamical symmetries. It addresses situations in which a prescribed symmetry is neither 
exact nor completely broken. When found, such intermediate symmetry structure can provide analytic 
solutions and quantum numbers for a portion of the spectra, thus offering considerable insight into 
complex dynamics. In other circumstances a PDS can serve as a convenient starting point for further 
improvements. 

As discussed in the present review, PDSs of various types have concrete applications to nuclear and 
molecular spectroscopy. Their empirical manifestations illustrate their ability to serve as a practical 
tool for calculation of observables in realistic systems. Hamiltonians with partial symmetries are not 
completely integrable nor fully chaotic. As such, they are relevant to the study of mixed systems with 
coexisting regularity and chaos, which are the most generic. Quantum phase transitions are driven by 
competing symmetries in the Hamiltonian. They provide a natural arena for PDSs, which incorporate 
such incompatible symmetries, and manage to survive at the critical points, in spite of the strong 
mixing. 

PDSs appear to be a common feature in algebraic descriptions of dynamical systems. They are 
not restricted to a specific model but can be applied to any quantal systems of interacting particles, 
bosons and fermions. The existence of PDS is closely related to the order of the interaction among 
the particles. The partial symmetry in question can be continuous (Lie-type), or discrete (point-group) 
and the associated dynamical algebra can involve a single or coupled algebraic structure. The examples 
of partial dynamical symmetries surveyed in the present review, involved purely bosonic or purely 
fermionic algebras. It is clearly desirable to extend the PDS notion to mixed systems of bosons and 
fermions, and explore the possible role of partial Bose-Fermi symmetries and partial supersymmetries. 

On phenomenological grounds, having at hand concrete algorithms for identifying and constructing 
Hamiltonians with PDS, is a valuable asset. It provides selection criteria for the a priori huge number 
of possible symmetry-breaking terms, accompanied by a rapid proliferation of free-parameters. This 
is particularly important in complicated environments when many degrees of freedom take part in the 
dynamics and upon inclusion of higher-order terms in the Hamiltonian. Futhermore, Hamiltonians with 
PDS break the dynamical symmetry (DS) but retain selected solvable eigenstates with good symmetry. 
The advantage of using interactions with a PDS is that they can be introduced, in a controlled manner, 
without destroying results previously obtained with a DS for a segment of the spectrum. These virtues 
greatly enhance the scope of applications of algebraic modeling of quantum many-body systems. 

On a more fundamental level, it is important to recall that dynamical systems often display simple 
patterns amidst a complicated environment. A representative example is the occurrence of both collec- 
tive and quasi-particle type of states in the same nucleus. It is natural to associate the "simple" states 
with a symmetry that protects their purity and special character. This symmetry, however, is shared 
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by only a subset of states, and is broken in the remaining eigenstates of the same Hamiltonian. It thus 
appears that realistic quantum many-body Hamiltonians can accommodate simultaneously eigenstates 
with different symmetry character. These are precisely the defining ingredients of a partial symmetry. 
In this context, PDS can offer a possible clue to the deep question of how simple features emerge from 
complicated dynamics. 

Underlying the PDS notion, is the recognition that it is possible for a non-invariant Hamiltonian 
to have selected eigenstates with good symmetry and good quantum numbers. In such a case, the 
symmetry in question is preserved in some states but is broken in the Hamiltonian (an opposite situation 
to that encountered in a spontaneously-broken symmetry). The PDS concept is, therefore, relevant to 
situations where selected states fulfill the predictions of a symmetry, which is otherwise known to be 
broken. Familiar examples of such a scenario include flavor symmetry in hadrons and chiral symmetry 
in nuclei. PDS may thus shed a new light on the related question of why, occasionally, symmetries seem 
to be obeyed beyond their domain of validity. 

The realistic attributes and fundamental aspects of partial symmetries, as portrayed in the present 
review, illuminate their potential useful role in dynamical systems and motivate their ongoing and 
future in-depth study. Much has been learnt but much more awaits to be explored and understood. 
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Appendix 



The normal-order form of the most general IBM Hamiltonian with one- and two-body interactions is 
given by 



H 



IBM 



e s s^s + e d d) ■ d + u (s^) 2 s 2 + u 2 s*(P ■ ds + v (s^) 2 d ■ d + H.c. 

+V 2 [fitrft . (JJ)(2) + HmC ] + C L ■ (dd) iL) . 

L=0,2,4 



H 



IBM 



e s N + u N(N - 1) + (e d - e s ) h d - (2u - u 2 + 2v ) (N - l)n d 
+ (u -u 2 + 2v + ^v 2 + \c Q + f c 2 + ^c 4 )n d (n d - 1) 



C\ 



0(5) 



4n„ 



2y/7 



t' 2 



;i8i) 



The corresponding multipole form, written in terms of the Casimir operators listed in Table [161, is given 
by 



+(w + 7^2 



C (6) - 5V 



+ + \ C 4 - \ C 2) 



a 



0(3) 



C SU (3) - ION 
- 6h, 



'182) 
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Table 16: Generators and Casimir operators, Cq, of algebras G in the IBM. Here h d = y/5U^, = 
yiOC/W, Q( 2 ) = n( 2 ) - ^UW, Q( 2 ) = n(2) + ^Z[/( 2 ). The operators U^, I^ 2 ), n„ are defined in 
Eq. ©. 
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